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I.  INTRODUCTION 

There  are  a  number  of  practical  applications  for  which  annular,  super¬ 
sonic  nozzles  are  used.  Among  these  are  turbofan  bypass  nozzles,  as  well 
as  unconventional  propulsion  nozzles  such  as  the  spike,  plug,  and  expansion- 
deflection  designs.  Tne  purpose  of  the  latter  group  of  nozzles,  which  all 
contain  centerbodies ,  is  to  obtain  improved  thrust  performance  over  the 
conventional  converging-diverging  configuration  at  off-design  operating 
conditions.  Another  application  of  annular  supersonic  nozzles  which  has 
recently  been  proposed  [1]*,  is  the  use  of  coaxial,  coflowing  supersonic 
streams  in  order  to  obtain  jet  noise  suppression.  The  advantage  of  this 
configuration  is  that  the  complex,  interacting  two-stream  flowfield  allows 
the  suppression  of  noise-generation  mechanisms  with  a  minimal  thrust 
penalty,  as  opposed  to  the  insertion  of  mechanical  devices  into  the  exhaust 
flow. 

However,  the  primary  motivation  for  this  work  is  the  supersonic- 
supersonic  ejector,  Fig.  1.1,  whereby  an  energetic,  coflowing,  supersonic, 
primary  stream  is  used  to  pump  a  low  pressure,  supersonic,  secondary  stream. 
This  device  has  applications,  for  example,  in  establishing  and  maintaining 
supersonic  flow  conditions  in  the  cavity  region  of  chemical  lasers  [2,3]  or 
as  a  means  of  augmenting  the  pressure  recovery  for  supersonic  wind  tunnel 
applications  [4-7]. 

From  previous  experience  with  the  subsonic-supersonic  ejector  [8,9],  for 
which  the  secondary  stream  is  subsonic,  it  is  known  that  the  overall  ejector 

lumbers  in  brackets  refer  to  entries  in  REFERENCES. 
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operating  characteristics  are  controlled  to  a  great  extent  by  the  interaction 
of  the  primary  and  secondary  streams  immediately  downstream  from  their 
point  of  confluence.  Thus,  any  detailed  investigation  of  the  supersonic- 
supersonic  ejector  flowfield  should  naturally  concentrate  on  this  "initial 
interaction  region."  For  example,. in  the  preliminary  investigation  of 
Guile,  et  al.  [10],  a  two-dimensional,  inviscid,  rotational  method  of 
characteristics  technique  was  used  to  analyze  the  flow  in  this  region  start¬ 
ing  at  the  ejector  inlet  station.  However,  in  this  adverse  pressure  gradient 
situation,  it  is  recognized  that  the  boundary  layer  and  viscous  shear  layer 
flows  near  the  primary-secondary  confluence  should  be  investigated.  The 
experiments  of  Mikkelsen,  et  al .  [2]  showed  severe  pressure  recovery  losses 
for  the  constant  area  supersonic-supersonic  ejector  due  to  separation  of  the 
secondary  stream  as  the  static  pressure  of  the  primary  was  increased  above 
that  of  the  secondary  at  their  confluence  point.  It  is  also  known  that  in 
the  hot-flow,  chemical  laser  application  the  boundary  layers  of  the  secondary, 
laser  flow  are  extremely  thick  at  the  ejector  station.  Therefore,  any 
detailed  calculations  of  the  primary-secondary  flowfield  interactions  should 
include  these  viscous  effects  and  preferably  should  be  initiated  upstream  of 
the  ejector  inlet  so  that  the  characteristics  of  the  boundary  layers  along 
the  walls  of  the  primary  and  secondary  streams  can  be  well  established 
before  the  confluence  point  is  reached. 

One  natural  place  to  start  these  calculations,  particularly  for  the 
primary  stream,  is  in  the  throat  region  of  the  supply  nozzles  using  an 
appropriate  transonic  throat  analysis.  Since  the  Reynolds  number  of  the 
primary  stream  is  generally  quite  large  and  since  the  streamwise  wall 
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pressure  gradient  in  the  throat  region  is  favorable,  the  inviscid  assumption 
may  be  made  in  the  development  of  the  throat  analyses.  The  results  can  then 
be  used  as  initial  conditions  to  start  the  coupled  invi scid/boundary  layer 
computations  which  proceed  downstream  through  the  remainder  of  the  ejector 
flowfi  el  d. 

For  both  the  chemical  laser  and  wind  tunnel  appl  icatio'ns  of  the 
supersonic-supersonic  ejector,  the  geometry  at  the  ejector  may  be  either 
two-dimensional  plane  or  axisymmetric.  Hence,  transonic  analyses  for 
nozzles  with  plane,  axisymmetric-no  centerbody,  or  annular  configurations 
may  be  required  to  initialize  the  ejector  flowfield  calculations.  In 
particular,  if  the  ejector  sketched  in  Fig.  1.1  is  taken  as  axisyrrmetric, 
the  primary  supply  nozzle  is  an  annular  one  with  a  small  relative  distance 
from  the  axis  of  symmetry  to  the  nozzle  when  non-dimensionalized  with 
respect  to  the  throat  height.  However,  if  the  primary  stream  is  somehow 
centrally  submerged  in  the  supersonic  secondary,  for  example  by  a  strut  sup¬ 
port  arrangement,  its  supply  nozzle  is  a  conventional  axisymmetric  one. 

Plane  configurations  for  the  supersonic-supersonic  ejector  have  also  been 
proposed  and  studied  [10]. 

The  objective  of  this  study  is  twofold.  First  it  is  desired  to  develop 
an  analytical  technique  which  provides  an  accurate  description  of  the 
transonic  throat  flowfield  for  a  wide  variety  of  two-dimensional  nozzle 
configurations  including  those  mentioned  previously.  Although  the  analysis 
is  to  be  applicable  to  as  wide  a  region  as  possible  in  the  nozzle  throat, 
its  main  intent  is  to  establish  an  initial  value  line  for  starting  the 
marching-type  computations,  e.g.,  method  of  characteristics,  for  the 
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supersonic  portion  of  the  flowfield.  Also,  since  the  implementation  of  such 
transonic  analyses  is  invariably  numerical,  it  is  desired  that  the  solution 
be  numerically  fast  and  reliable  so  that  parametric  studies  or  iterative 
calculations,  as  might  occur  in  nozzle  design  situations,  can  be  made  in  an 
efficient  and  low-cost  manner.  The  second  objective  is  to  provide  further 
insight  into  the  characteristics  of  transonic  throat  flows  by  an  experi¬ 
mental  investigation  of  a  number  of  configurations  including  conventional 
axisynmetric  and  annular  ones.  Comparison  of  the  theoretical  and  experi¬ 
mental  results  then  allows  conclusions  to  be  drawn  regarding  the  limits  of 
applicability  of  the  analysis  and  suggests  possible  refinements  of  it. 

As  will  be  discussed  further  in  the  next  chapter,  no  analytical  tech¬ 
niques  are  known  to  exist  which  contain  all  of  the  desired  features  just 
outlined.  In  addition,  experimental  information  concerning  transonic 
throat  flows  is  extremely  sparse,  particularly  for  axisymmetric  geometries. 
This  investigation  is  intended  to  provide  additional  information  in  these 
areas . 


ej ec tor- di f fuser  con  figuration 


II.  LITERATURE  REVIEW 


The  subject  of  transonic  flow  in  nozzles  and  ducts  has  been  studied 
for  well  over  a  hundred  years,  dating  from  the  converging  nozzle  studies  of 
Navier  [11]  in  1829  and  St.  Venant  and  Wantzel  [12]  in  1839.  During  this 
period,  many  techniques  have  been  developed  for  investigating  nozzle  threat 
flows  as  the  state  of  understanding  of  this  topic  has  been  advanced.  The 
large  majority  of  the  contributions  in  this  area  are  theoretical  in  nature 
and  apply  to  the  conventional  two-dimensional  plane  and  axisymmetric  con¬ 
figurations,  since  study  of  annular  nozzles  has  been  motivated  by  the 
relatively  recent  development  of  devices  such  as  the  plug  and  expansion- 
deflection  nozzles,  the  supersonic-supersonic  ejector,  etc.  Because  the 
annular  configuration  is  of  prime  importance  in  this  investigation,  the 
purpose  of  the  present  chapter  is  to  review  those  papers  dealing  explicitly 
with  this  configuration  and  to  summarize  the  major  contributions  to  the 
general  methods  which  have  been  developed  for  throat  flow  analysis.  The 
survey  of  Hall  and  Sutton  [13]  provides  a  useful  review  of  early  work  while 
the  more  recent  surveys  of  Flack  and  Thompson  [14]  and  Brown  and  Hamilton 
[15]  consider  newer  developments  in  this  area. 

A.  ANALYTICAL 

One  of  the  fundamental  difficulties  involved  with  the  theoretical  solu¬ 
tion  of  most  fluid  flow  problems  is  the  basic  nonlinear  nature  of  the  govern¬ 
ing  equations.  Further  complicating  the  analysis  of  steady  transonic  flows 
is  the  fact  that  in  the  subsonic  portions  of  the  flowfield  the  equations  are 
of  elliptic  type  while  in  the  supersonic  portions  they  are  of  hyperbolic 
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type.  The  mixed,  nonlinear  nature  of  the  equations  describing  transonic 
nozzle  throat  flows  makes  the  solution  of  these  problems  by  both  analytical 
and  numerical  techniques  difficult  at  best.  However,  several  methods  of 
varying  degrees  of  accuracy  and  usability  have  been  developed.  Those  to 
be  reviewed  in  this  section  have  been  divided  into  the  following  categories: 
indirect  techniques,  series  expansion  methods,  time  dependent  numerical 
techniques,  and  other  methods. 

1.  Indirect  Techniques 

In  the  indirect  or  inverse  methods,  the  velocity  distribution  is 
assumed  along  a  reference  streamline,  such  as  the  nozzle  centerline,  the 
flowfield  is  calculated  by  some,  usually  approximate,  technique,  and  then 
in  the  inviscid  approximation,  any  other  streamline  can  be  taken  as  a 
nozzle  wall.  Therefore,  indirect  techniques  can  be  applied  only  to  the 
nozzle  design  problem  and  cannot,  without  iteration,  be  applied  to  the 
direct  problem  of  calculating  the  flowfield  in  a  nozzle  with  given  wall 
contours. 

Among  the  earliest  attempts  at  analyzing  the  flow  in  the  throat  of  a 
plane  converging-diverging  nozzle  was  that  of  Meyer  [16].  He  assumed  that 
the  velocity  along  the  axis  increased  linearly  and  substituted  a  double 
power  series  expansion  in  the  coordinates  into  the  perturbation  potential 
equation.  The  coefficients  up  to  and  including  the  sixth  order  terms  were 
obtained  from  which  the  characteristics  of  the  flowfield  were  deduced. 

Another  early  solution  technique  which  was  applied  to  nozzle  flow 
analyses  is  the  hodograph  method  developed  successively  by  Lighthill  [17,18], 
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Frank!  [19],  and  Cherry  [20-22],  In  this  method  the  hodograph  equations, 
which  employ  the  velocity  components  as  independent  variables,  are  used. 
These  equations  are  linear,  but  the  solution  is  complicated  by  the  fact 
that  over  a  portion  of  the  supersonic  flowfield  the  mapping  from  the  hodo¬ 
graph  plane  to  the  physical  plane  is  not  one-to-one.  Cherry  [20]  was  able 
to  devise  a  transformation  which  provides  a  si ngl e-val ued  solution  over  the 
entire  plane  and  later  [21,22]  showed  how  superposition  could  be  used  to 
obtain  solutions  for  realistic  nozzle  flows.  Major  drawbacks  of  the  hodo¬ 
graph  method  are  that  it  is  applicable  only  to  the  planar  case  and  that 
practical  calculations  appear  to  be  quite  complicated. 

An  early  publication  dealing  with  annular  nozzles  is  that  of  Lord 
[23]  who  considered  a  supersonic  nozzle  with  a  cylindrical  centerbody.  In 
his  inverse  technique,  two  orders  of  solution  were  obtained  by  assuming  that 
the  velocity  along  the  cylinder  surface  increases  linearly  in  the  streamwise 
direction  for  the  first  order  solution  and  quadratical ly  for  the  second 
order.  As  noted  by  Lord,  the  solutions  for  the  annular  case  are  much  more 
complicated  than  for  the  limiting  cases  of  axisymmetric  and  plane  nozzles 
for  which  the  centerbody  diameter  approaches  zero  and  the  throat  diameter, 
respectively.  Lord  used  his  throat  flowfield  results  to  obtain  a  starting 
line  for  method  of  characteristics  nozzle  design  computations. 

One  of  the  more  well  known  indirect  techniques  was  developed  by  Hopkins 
and  Hill  [24]  who  used  the  Friedrichs  formulation  to  obtain  the  flowfields 
in  axisymmetric  nozzles  with  a  small  wall  radius  of  curvature  at  the  throat. 
The  dependent  variables  were  expanded  as  power  series  in  the  stream  function 
so  that  the  flowfield  could  be  found  from  the  prescribed  velocity 
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distribution  along  the  axis  by  solving  for  the  series  coefficients.  A  corre¬ 
lation  was  developed  between  a  one-dimensional  "reference  boundary"  and  the 
computed  boundary  to  facilitate  the  solution  of  the  direct  problem  of 
analyzing  nozzles  with  specified  wall  contours.  The  authors  found  that  the 
most  significant  geometric  parameter  affecting  the  flow  in  the  throat  is  the 
wall  radius  of  curvature  at  the  throat  and  that  the  inlet  angle  of  the  con¬ 
vergent  section  is  not  important  unless  this  radius  of  curvature  is  less 
than  the  throat  radius.  In  a  later  publication  [25],  Hopkins  and  Hill 
extended  their  method,  with  corrections  by  Morden  and  Farquhar  [26],  to  the 
case  of  unconventional,  annular  nozzles  inclined  at  arbitrary  angles  to  the 
axis  of  symmetry.  Their  method  assumes,  however,  the  existence  of  a 
straight  reference  streamline  along  which  the  velocity  is  specified;  this 
is  not  necessarily  the  case  in  the  direct  problem.  Also,  because  of  the 
asymmetry  of  the  streamlines  about  the  reference  streamline,  no  correlations 
were  obtained  between  reference  and  computed  boundaries.  For  analysis  of 
given  nozzles,  iterations  on  a  computer  are  suggested. 

Norton  [27]  developed  a  somewhat  different  inverse  technique  for  con¬ 
ventional  axisymmetric  nozzles.  In  his  method,  the  Euler  equations,  includ¬ 
ing  the  effects  of  rotational ity  and  swirl,  are  first  transformed  by  intro¬ 
ducing  the  stream  function  and  a  stretched  axial  coordinate;  these 
transformed  equations  are  then  numerically  integrated  outward  from  the  given 
velocity  distribution  along  the  axis.  Norton  presents  some  calculated 
results,  which  required  approximately  two  minutes  of  Univac  1108  time,  but 
emphasizes  that  the  method  is  not  suited  to  the  problem  of  analyzing  given 
nozzle  contours. 


In  a  novel  application  of  the  indirect  method,  Van  Tuyl  [28]  has 
obtained  solutions  for  both  axisymmetric  and  planar  nozzles.  The  stream 
function  and  density  are  expanded  as  Taylor  series  in  the  neighborhood  of  a 
point  on  the  nozzle  axis  along  which  the  velocity  distribution  is  specified. 
The  coefficients  of  these  series  are  found  by  substitution  into  the  govern¬ 
ing  irrotational i ty  and  Bernoulli  equations,  and  series  for  other  quantities 
of  interest  are  also  determined.  Pad£  fractions  are  then  used  to  obtain 
convergence  when  the  resulting  power  series  diverge  and  to  accelerate  con¬ 
vergence  when  they  converge.  Van  Tuyl  carried  out  sample  calculations  for 
a  specific  axisymmetric  nozzle  using  a  digital  computer  both  to  obtain  the 
terms  of  the  series  and  to  form  the  Pad£  approximants .  The  computations 
took  21  seconds  on  a  CDC  6400  when  25  terms  in  the  power  series  were  retained. 

An  inverse  scheme  for  analyzing  the  flow  in  conventional  plane  and  axi¬ 
symmetric  nozzles  has  been  reported  by  Klopfer  and  Holt  [29].  In  their 
formulation,  the  von  Mises  transformation  is  applied  to  the  continuity  and 
Euler  equations  so  that  the  dependent  variables  are  taken  as  the  pressure, 
the  transverse  coordinate,  and  the  ratio  of  the  transverse  to  streamwise 
velocity  components.  The  resulting  equations  are  numerically  integrated 
using  the  prescribed  pressure  distribution  along  the  nozzle  axis.  Results 
for  both  plane  inverse  and  axisymmetric  direct  problems  are  presented,  and 
good  agreement  is  obtained  with  the  data  of  Cuffel ,  et.  al.  [30]  in  the 
axisymmetric  case.  The  direct  problem  is  solved  iteratively  as  a  sequence 
of  inverse  problems. 

Ishii  [31]  has  recently  proposed  an  inverse  method  for  nozzle  design. 

In  this  technique,  the  governing  equations  for  inviscid  flow  are  transformed 
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so  as  to  use  the  streamline  and  orthogonal  trajectories  as  independent  vari¬ 
ables.  The  flow-field  is  divided  into  streamtubes  and  integration  of  the 
equations  is  carried  out  numerically  starting  from  the  prescribed  velocity 
distribution  along  the  nozzle  axis. 

As  should  be  clear  from  the  preceding  discussion,  the  major  disadvan¬ 
tage  of  the  indirect  methods  is  their  inability  to  deal  with  the  problem  of 
analyzing  the  flow  in  a  nozzle  of  a  given  configuration.  Only  Hopkins  and 
Hill  [24]  have  devised  a  scheme  for  avoiding  the  iterations  which  are  inevi¬ 
table  when  applying  an  inverse  technique  to  the  direct  problem.  As  discussed 
by  Brown  and  Hamilton  [15]  and  Klopfer  and  Holt  [29],  another  objection  to 
these  methods  is  the  potentially  poorly  posed  nature  of  the  problem  resulting 
from  the  specification  of  the  velocity  along  the  centerline  in  the  subsonic 
region.  When  the  governing  equations  are  solved  numerically,  this  specifi¬ 
cation  is  known  to  have  a  destabilizing  effect  on  the  computations  (Hadamard 
instability). 

2.  Series  Expansion  Methods 

The  series  expansion  techniques  have  proven  to  be  very  popular  for 
obtaining  the  flowfields  in  the  throats  of  supersonic  nozzles.  In  these 
methods,  the  dependent  variables  are  expanded  in  series,  and  the  unknown 
constants  and  functions  are  then  determined  by  substitution  into  the  govern¬ 
ing  equations  and  boundary  conditions.  These  techniques  can  generally  be 
considered  to  provide  approximate  solutions  for  the  problem  of  analyzing  the 
flow  in  a  nozzle  of  a  given  configuration. 

One  of  the  first  attempts  at  solving  the  direct  problem  of  flow  in  a 
planar  nozzle  was  made  by  Taylor  [32]  in  1930.  He  expanded  the  compressible 
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perturbation  potential  in  a  double  power  series  in  the  coordinates  and  sub¬ 
stituted  it  into  the  corresponding  potential  equation.  By  using  the 
boundary  condition  that  the  nozzle  wall  is  a  streamline  in  the  inviscid 
sense,  the  coefficients  in  the  series  were  obtained.  Taylor  calculated  an 
example  for  a  nozzle  with  circular  walls  and  a  dimensionless  radius  of 
curvature1 ,  R  ,  of  4.  , 
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Somewhat  later,  Sauer  [33]  presented  series  solutions  for  the  flow  in 
both  axisymmetric  and  planar  supersonic  nozzles.  In  this  well  known  tech¬ 
nique,  the  compressible  potential  equation  was  first  simplified  by  retaining 
only  the  dominant  terms,  and  the  perturbation  potential  function  was  then 
expanded  as  a  power  series  in  the  transverse  coordinate,  y,  with  unknown 
coefficient  functions  of  the  streamwise  coordinate,  x.  By  substituting  the 
expansion  into  the  simplified  potential  equation  and  relating  the  resulting 
flowfield  to  the  curvature  of  the  nozzle  walls,  the  unknown  functions  and 
constants  were  obtained.  However,  since  the  series  was  terminated  after  the 
first  two  terms,  the  results  are  applicable  only  to  nozzles  with  small  wall 
curvature,  i.e.,  large  R^  . 

At  about  the  same  time,  a  somewhat  different  solution  technique  for  the 
conventional  configurations  was  proposed  by  Oswatitsch  and  Rothstein  [34]. 
Rather  than  obtaining  solutions  for  nozzles  whose  centerline  velocity  dis¬ 
tribution  is  approximated  to  be  linear  as  in  Sauer's  method,  any  centerline 
distribution  is,  in  principle,  possible  in  Oswatitsch  and  Rothstein's  method. 

throughout  this  investigation,  the  dimensionless  radius  of  curvature,  R  , 
is  defined  as  the  throat  wall  radius  of  curvature  divided  by  the  throatc 
half-height.  For  conventional  axisymmetric  nozzles  the  half-height  is 
taken  as  the  throat  radius,  while  for  annular  nozzles  it  is  taken  as  the 
distance  between  the  inner  and  outer  throat  wall  locations  in  the 
meridional  plane. 


however,  since  it  is  determined  by  the  iterative  solution  of  a  differential 
equation  which  results  from  the  continuity  equation  for  a  given  nozzle  con¬ 
tour,  this  technique  is  much  more  difficult  to  use  than  Sauer's  closed  form 
solutions.  In  addition,  Oswatitsch  and  Rothstein's  method  gives  anomalous 
results  for  nozzles  which  have  a  small  wall  radius  of  curvature. 

Hall  [35]  also  carried  out  an  expansion  solution  for  the  flow  in  the 
throat  of  axi symmetri c  and  plane  nozzles.  In  this  investigation,  a  careful 
order  of  magnitude  analysis  of  the  various  dependent  and  independent  vari¬ 
ables  was  made,  and  the  transonic  perturbation  velocity  components  were 
expanded  in  appropriate  series  using  .  =  1  as  the  expansion  parameter. 

By  substituting  these  series  into  the  wall  boundary  conditions  and  the 
governing  irrotational ity  and  gas  dynamic  equations,  the  formulations  for 
the  various  orders  of  solution  in  the  expansion  technique  were  obtained. 
Closed  form  expressions  for  the  first  three  orders  were  determined  by  using 
solution  forms  suggested  by  the  boundary  conditions.  The  first  order  solu¬ 
tions  are  identical  to  those  of  Sauer  [33],  and  Hall  also  found  that  to  the 
first  three  orders  the  approximate  solutions  for  nozzles  with  circular, 
parabolic,  and  hyperbolic  arc  wall  contours  coincide.  Nozzles  for  which 
the  wall  radius  of  curvature  is  less  than  the  throat  half-height,  Re<l, 
cannot  be  analyzed  with  this  technique  since  this  results  in  expansion 
parameters  greater  than  unity.  In  fact  the  smallest  value  of  Rc  used  by 
Hall  in  his  example  calculations  is  10/3. 

This  same  solution  technique  was  later  extended  by  Moore  and  Hall  [36] 
to  the  case  of  annular  nozzles  inclined  arbitrarily  to  the  axis  of  symmetry. 
However,  as  pointed  out  by  Thompson  and  Flack  [37],  the  coordinate  system 


I 

I 


I. 


In  a  con- 


employed  in  this  study  is  an  overspecified,  inconsistent  one. 
panion  report,  Moore  [381  used  •  =  R  ‘‘as  the  expansion  parameter  to 
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construct  approximate  solutions  for  the  flow  in  asymmetric  planar  nozzles. 
The  first  three  terms  of  the  series  solution  were  presented  along  with  the 
results  for  two  special  cases,  namely  the  conf igurations  in  which  the  two 
nozzle  walls  have  unequal  curvature  and  in  which  the  nozzle  wall  contour 
is  asymmetric  about  the  throat  plane. 

In  a  series  of  articles  [39-42],  Szaniawski  has  investigated  transonic 
flow  in  nozzle  throats.  In  [39]  and  [40],  the  governing  equations  are 
derived  and  the  existence  and  nature  of  solutions  for  the  flow  of  a  viscous 
heat  conducting  gas  in  the  throat  region  of  two-dimensional  plane  nozzles 
are  investigated.  The  two  later  papers  [41,42]  are  concerned  with  finding 
solutions  under  the  usual  assumptions  of  inviscid,  irrotational  flow.  A 
solution  is  developed  for  a  particular  class  of  planar  nozzle  wall  contours 
in  [41].  The  potential  function  is  expanded  as  a  power  series  in  the  trans 
verse  coordinate  with  undetermined  coefficient  functions  of  the  streamwise 
coordinate.  Recurrence  relations  for  these  functions  are  then  determined 
by  substitution  into  the  potential  equation  and  the  boundary  conditions.  A 
similar  method  is  used  in  [42]  to  determine  solutions  for  both  plane  and 
axisymmetric  conventional  nozzles,  but  an  auxiliary  system  of  curvilinear 
coordinates  is  employed.  Both  solution  techniques  are  applicable  only  to 
nozzles  which  have  a  large  throat  radius  of  curvature. 

In  a  somewhat  different  but  closely  related  area,  Sichel  [43,44]  has 
studied  the  flow  in  two-dimensional  nozzle  throats  at  operating  conditions 
near  the  choking  point.  The  transition  from  symmetric  Taylor  flow  to 
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continuously  accelerating  Meyer  flow  was  investigated  by  including  the 
effect  of  longitudinal  viscosity  in  the  governing  equations.  Based  on  this 
formulation,  a  similarity  solution  was  found  which  provides  a  smooth  tran¬ 
sition  between  the  regimes  and  which  shows  the  initial  stages  of  shock 
formation. 

kliegel  and  Quan  [45]  used  a  method  similar  to  Hall's  for  predicting 
plane  and  axi symmetri c  nozzle  flows  except  that  a  different  stretched  axial 
coordinate  was  employed.  The  velocity  components  were  expanded  as  series 
in  •  =  R  1 ,  and  two  orders  of  corrections  to  the  one-dimensional  solution 
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were  obtained.  The  authors  conclude  that  the  resulting  solution,  which  is 
claimed  to  be  valid  throughout  the  entire  nozzle,  should  be  limited  to 
nozzles  which  have  a  large  wall  radius  of  curvature,  Rc_>2,  and  that  Hall's 
transonic  solution  can  be  obtained  from  it  through  a  coordinate  transforma¬ 
tion  and  reordering  of  the  terms.  The  cases  of  unchoked  nozzle  flows  and 
multistream  nozzle  expansions,  as  would  occur  in  film  cooled  rocket  engines, 
are  also  discussed. 

The  expansion  solutions  presented  so  far  are  all  limited  to  nozzles 
with  gentle  wall  curvatures.  In  an  attempt  to  analyze  the  case  of  small 
wall  radius  of  curvature  for  axisymmetric  nozzles,  Kliegel  and  Levine  [46] 
presented  a  solution  utilizing  e  =  IR  +1 )"'  as  the  expansion  parameter. 

c 

Since  this  parameter  is  less  than  unity  for  all  values  of  the  radius  of 
curvature,  the  resulting  series  should  have  superior  convergence  properties 
in  the  limit  of  small  R  .  The  authors  contend  that  this  solution  is  the 
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one  obtained  when  the  solution  in  toroidal  coordinates,  which  are  convenient 
for  circular  arc  wall  contours,  are  transformed  back  to  cylindrical 


coordinates.  However,  in  a  later  publication  Levine  and  Coats  [47]  conclude 
that  this  contention  is  false  and  also  that  'the  pmnnsed  series  do  not 
satisfy  the  differential  equations  of  motion  in  cylindrical  coordinates . 

In  reality,  the  "solution"  proposed  in  { 46 j  is  simply  Hall's  solution  in 
•  =  R  1  transformed  to  a  series  in  •  =  (R  +1 )  1  such  that  the  two  are 
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equivalent  in  the  limit  of  large  radius  of  curvature,  R  .  xlieqel  and 
Levine  also  corrected  some  errors  which  appeared  in  Hall's  third  order 
solution  for  axi symmetri c  nozzles,  although  the  present  author  has  found 
their  correction  of  the  discharge  coefficient  to  be  in  error. 

The  conclusions  of  Levine  and  Coats  [47]  just  mentioned  were  the 
result  of  an  extensive  investigation  of  axisymmetric  nozzle  throat  flows. 

The  transonic  equations  of  motion,  including  the  effect  of  variable 
specific  heat  ratio,  were  solved  in  toroidal  coordinates  using  a  combina¬ 
tion  of  a  small  parameter  expansion  and  a  double  series  in  the  coordinates. 
Although  a  number  of  expansion  parameters  was  tried,  none  of  the  resulting 
solutions  was  found  to  be  convergent  for  small  values  of  Rc  .  The  authors 
hypothesized  that  this  is  because  the  expansion  solutions  assume  that  the 
local  wall  geometry,  in  particular  the  throat  radius  of  curvature,  completely 
determines  the  flowfield  in  the  throat  and  that  this  premise  breaks  down  as 
R  ->0.  In  addition,  the  effect  of  variable  specific  heat  ratio  was  found  to 
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be  negligible  in  the  transonic  region. 

Taulbee  and  Boraas  [483  have  developed  a  Hall-type  series  solution  for 
transonic  flow  in  the  throat  region  of  axisymmetric  nozzles  with  a  nonuni¬ 
form  stagnation  temperature  distribution  at  the  inlet.  The  governing  con¬ 
tinuity  and  Euler  equations  were  first  transformed  using  the  stream  function 


as  an  independent  variable,  and  expansions  for  the  dependent  variables  were 
then  assumed  using  an  order  of  magnitude  analysis.  Two  orders  of  solutions 
were  obtained,  and  the  results  indicate  that  the  sonic  line  location  is 
essentially  unaffected  by  the  nonuni formi ty  of  the  stagnation  temperature. 
However,  as  demonstrated  by  the  calculated  example  presented  in  [49],  the 
flowfield  details  such  as  the  flow  angle  sensitive  to  the  total  temper¬ 
ature  di strioution.  In  a  later  paper,  Boraas  [50]  included  the  effects  of 
nonuni  form  total  temperature  and  composition  in  an  analysis  of  transonic 
throat  flow  in  axisymmetric  nozzles.  First  and  second  order  solutions  were 
again  found  (see  Ishii  [51]  for  corrected  equations),  and  it  is  stated  that 
in  this  case  the  sonic  line  location  does  depend  on  the  nonuniformities  of 
the  properties  upstream  of  the  throat  although  no  examples  are  presented. 

In  both  of  the  rotational  analyses  just  cited  the  expansion  parameter  is 
taken  as  e  =  R^'  so  that  the  solutions  are  not  applicable  to  nozzles  with 
small  radius  of  curvature  of  the  throat  wall. 

Apparently  unaware  of  the  earlier  work  of  lord  [23],  Smithey  and  Naber 
[52]  presented  a  solution  for  the  sonic  line  in  an  axisymmetric  nozzle  with 
a  concentric  inner  cylinder.  The  solution  was  obtained  by  hypothesizing  a 
trial  form  suggested  by  the  Sauer  [33]  solutions  for  the  limiting  cases  of 
axisynmetric  and  plane  nozzles.  Substitution  into  the  transonic  gas 
dynamic  equation  and  use  of  the  boundary  conditions  then  allowed  determina¬ 
tion  of  the  unknown  functions  and  constants.  However,  since  the  resulting 
solution  satisfies  the  irrotational ity  condition  only  for  the  limiting  con¬ 
figurations,  its  applicability  to  the  annular  case  of  interest  is 
questionable. 


Thompson  and  Flack  [14,37,53]  have  performed  a  very  extensive  investi¬ 
gation  of  the  flow  in  the  throat  region  of  asymmetric  plane  and  annular 
ax i symmetric  supersonic  nozzles.  In  their  theoretical  analysis,  the  method 
of  Hall  was  used  with  an  expansion  parameter  which  reduces  to  •  =  (R  +•  )"' 
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where  R  is  an  average  dimensionless  radius  of  curvature  for  the  two  walls 
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and  r  may  be  taken  as  zero  or  one.  Three  solution  orders  were  obtained, 
and  the  authors  noted  that  significant  improvement  in  the  convergence 
properties  of  the  solution  for  small  radii  of  curvature  resulted  from  using 
r  =  1.  However,  because  of  the  manner  in  which  the  axisymmetric  term  in 
the  gas  dynamic  equation  was  expanded,  the  annular  solutions  can  be  applied 
only  to  configurations  for  which  the  distance  from  the  axis  of  symmetry  to 
the  nozzle  throat  is  large  when  non-dimensional ized  with  respect  to  the 
throat  half-height,  i.e.,  annular  nozzles  which  approach  the  plane  configur¬ 
ation.  Their  solution  cannot  be  used  for  conventional  axisymmetric  nozzles, 
and,  in  addition,  it  appears  that  it  gives  planar  results  for  annular 
nozzles  which  are  not  inclined  with  respect  to  the  axis  of  symmetry. 
Experiments  utilizing  both  static  wall  pressure  and  laser  doppler  veloci- 
meter  measurements  were  performed  for  a  series  of  plane  nozzle  throat  con¬ 
figurations;  good  agreement  was  found  between  the  analytically  and  experi¬ 
mentally  determined  sonic  lines.  These  experiments  will  be  discussed  in 
more  detail  in  a  later  section. 

Recently,  Ishii  has  published  the  results  of  analytical  investigations 
of  the  flow  in  axisynmetri c  nozzle  throats  with  rotational ity  effects 
included.  In  [54]  throat  flowfields  with  rate  processes  and  nonuniform 
property  distributions  are  considered;  see  also  the  comment  by  Boraas  [55] 


and  the  rebuttal  by  Ishii  {51].  The  expansion  method  of  Hall  was  used  with 
the  parameter  taken  as  e  =  Rc  1 ,  and  one  term  of  the  series  is  presented 
for  frozen  and  equilibrium,  uniform  flows.  The  transonic  flow  of  a  vibra- 
tionally  relaxing  or  chemically  reacting  gas  in  a  converging-diverging 
nozzle  is  investigated  in  [56],  Again  the  Hall  method  was  employed,  but  in 
this  study  the  expansion  parameter  was  chosen  as  c  =  (R^+l)"1. 

The  main  disadvantage  of  the  series  expansion  techniques  is  that  they 
can  generally  be  applied  only  to  a  given  class  of  configurations.  In 
particular,  most  of  the  analyses  mentioned  above  have  been  developed  only 
for  the  conventional  two-dimensional  plane  and  axisymmetri c  configurations. 
Also,  for  all  of  these  analyses,  there  is  a  lower  limit  on  the  throat  wall 
radius  of  curvature  below  which  the  analyses  give  unreasonable  results.  On 
the  other  hand,  since  the  series  methods  generally  lead  to  closed  form 
expressions,  the  implementation  of  these  solutions  is  quite  straightforward. 

3.  Time  Dependent  Numerical  Techniques 

In  order  to  avoid  the  mixed  nature  of  the  governing  equations  for 
steady  transonic  flow,  the  time  dependent  methods  utilize  the  unsteady  form 
of  the  equations  which  are  hyperbolic  for  all  speed  regimes.  The  resulting 
initial  value  problem  is  solved  numerically  by  continuing  the  computations 
in  the  time  domain  from  an  arbitrarily  assumed  set  of  initial  conditions 
until  the  steady  state  is  asymptotically  approached. 

Saunders  {571  apparently  was  the  first  to  apply  this  technique  to 
internal  transonic  flows.  In  this  study  the  two-step  Lax-Wendroff  differenc¬ 
ing  scheme  was  employed  to  integrate  the  equations  for  flow  in  axisymmetric 
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nozzles  with  uniform  inlet  conditions  assumed.  Good  agreement  was  found 
between  the  numerical  results  and  the  wall  static  pressure  data  of  Back, 
et  al.  [58]  for  a  converging-diverging  nozzle.  The  computation  time  was 
reported  as  approximately  45  minutes  on  a  CDC  3200  computer. 

Migdal,  Klein,  and  Moretti  [59]  utilized  a  nonconservative  form  of  the 
basic  equations  together  with  a  stretching  of  the  axial  coordinate  to 
obtain  the  throat  flowfield  in  axisymmetric  Laval  nozzles.  An  example 
calculation  was  performed  for  a  configuration  with  a  dimensionless  throat 
wall  radius  of  curvature  of  0.625  and  inlet  and  exit  angles  of  45°  and  15°, 
respectively.  The  agreement  between  these  results  and  the  wall  pressure 
data  of  Back,  et  al .  [58]  is  excellent.  The  running  time  for  this  case 
was  less  than  five  minutes  on  an  IBM  360/75. 

By  using  the  two-step  Law-Wendroff  scheme  developed  by  Saunders  [57], 
Wehofer  and  Huger  [60]  investigated  the  effect  of  nonuniform  inlet  condi¬ 
tions  on  axisymmetric  convergent-divergent  nozzle  flows.  These  authors 
also  considered  the  case  of  an  axisymmetric  converging  nozzle  with  a  coni¬ 
cal  inlet  section  including  the  influence  of  the  operating  pressure  ratio 
and  the  associated  free  jet  structure.  The  computations  were  very  time 
consuming,  however,  taking  from  two  to  five  hours,  depending  on  the  mesh 
size,  on  an  IBM  360/50. 

Brunei  1  [61]  used  a  somewhat  different  technique  to  study  the  transonic 
flowfields  in  planar,  supersonic  nozzles.  In  his  method,  a  modified  upwind 
differencing  scheme  was  utilized  to  integrate  a  conservative  form  of  the 
governing  equations.  The  results  appear  to  be  reasonable  for  operating  con¬ 
ditions  in  which  the  nozzle  is  not  choked,  but  problems  were  encountered 
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for  the  choked  conditions  which  are  of  primary  interest.  In  the  latter 
case,  the  constant  Mach  number  lines  contain  slope  discontinuities  which 
Brunei  1  hypothesizes  are  due  to  numerical  pressure  disturbances  that  are 
propagated  from  the  supersonic  wall  region  into  the  subsonic  region. 
Experiments  were  also  conducted  using  an  annular  nozzle  with  a  cylindrical 
centerbody  and  small  throat  gap,  but  in  all  cases  only  unchoked  operating 
conditions  were  studied. 

In  order  to  predict  the  flowfields,  and  especially  the  discharge 
coefficients,  for  axisymmetric  and  plane  nozzles  with  small  throat  radius 
of  curvature,  Laval  [62,631  used  a  time  dependent  method  similar  to  that  of 
Saunders  [57] .  In  this  treatment,  however,  an  axial  coordinate  stretching 
was  employed  together  with  the  introduction  of  a  pseudo-viscosity  term  to 
avoid  stability  problems.  Calculations  were  carried  out  for  nozzles  with 
radius  of  curvature  values  down  to  0.1,  and  very  good  agreement  between 
computed  and  measured  values  of  the  flow  coefficient  was  noted.  Laval 
states  that  the  method  can  be  applied  to  annular  configurations,  although 
no  such  examples  are  presented  in  the  references  cited  above.  The  calcula¬ 
tions  required  approximately  two  hours  of  IBM  360/50  time. 

Serra  [64]  has  also  developed  a  transient  technique  for  the  determina¬ 
tion  of  nozzle  flows.  In  his  method,  the  one-step  Law-Wendroff  difference 
scheme  is  incorporated  with  a  numerical  damping  term  so  as  to  stabilize  the 
numerical  procedure  in  the  vicinity  of  shock  waves.  Also,  since  mathemati¬ 
cal  overspecification  of  the  entrance  flow  conditions  caused  accuracy  diffi¬ 
culties  in  earlier  numerical  studies,  the  problem  was  formulated  to  establish 
realistic  conditions  at  the  entrance.  The  computed  results  compare  well 
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with  experimental  data  and  previous  theories  for  axisyninetri c  and  planar 
nozzles.  An  example  is  also  included  of  the  flow  in  an  annular  turbofan 
passage.  The  computation  time  on  a  Univac  1108  computer  was  80  minutes 
for  an  axisymmetric  nozzle  problem. 

In  an  effort  to  reduce  the  long  computation  times  associated  with  time- 
dependent  methods,  Cline  [65-68]  has  developed  an  efficient  program  for 
calculating  a  wide  variety  of  compressible,  internal  flows.  The  second- 
order  accurate,  MacCormack  finite  difference  method  is  utilized  to  solve  the 
governing  equations  which  are  left  in  a  nonconservati ve  form.  A  second- 
order,  reference  plane  characteristic  scheme  is  used  for  the  boundary 
points,  and  an  explicit  artificial  viscosity  is  included  for  shock  calcula¬ 
tions.  The  most  recent  version  [68]  of  the  code  solves  the  Navier-Stokes 
rather  than  the  Euler  equations  so  that  viscous  flows  with  Reynolds  numbers 
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up  to  about  10  can  be  analyzed.  Several  examples  are  presented  including 
inviscid  converging,  converging-diverging,  and  annular  nozzle  flows  as  well 
as  a  viscous  converging-diverging  nozzle  calculation.  Good  agreement 
between  the  computations  and  experimental  data  is  obtained  with  computation 
times  ranging  from  29  seconds  on  a  CDC  6600  to  7  minutes  on  a  CDC  7600  for 
the  example  problems.  These  times  are  obtained,  however,  at  the  expense  of 
using  coarse  mesh  spacing,  particularly  in  the  inviscid  cases. 

The  major  advantage  of  the  time-dependent  numerical  techniques  is  that 
if  the  codes  are  developed  in  a  sufficiently  general  manner  they  may  be 
employed  in  the  analysis  of  a  variety  of  transonic  internal  flows.  A 
serious  drawback  to  their  use,  however,  is  the  long  computation  times 
required.  This  means  that  parametric  studies  or  iterative  calculations,  such 
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as  might  be  necessitated  in  a  design  procedure,  would  be  numerically  time 
consuming  and  expensive. 

4.  Other  Methods 

A  number  of  other  methods  have  also  been  used  for  nozzle  flow  analysis. 
One  such  technique  is  the  method  of  integral  relations  developed  originally 
by  Dorodnitsyn  [69].  In  the  application  of  this  method,  the  flowfield  is 
divided  into  streamwise  strips  across  which  the  governing  partial  dif¬ 
ferential  equations  are  integrated  using  assumed  profiles  for  the  dependent 
variables.  The  resulting  ordinary  differential  equations  are  then  solved 
to  obtain  the  coefficients.  Holt  [70]  formulated  the  one  strip  and  two 
strip  approximations  for  both  the  inverse  and  direct  problems  of  flow  in 
planar  nozzles.  Results  of  a  one  strip,  inverse  calculation  are  presented 
and  are  found  to  agree  well  with  the  hodograph  results  of  Cherry  [21]  in  the 
supersonic  region  but  to  diverge  seriously  from  them  in  the  subsonic  region. 
Liddle  and  Archer  [71]  carried  out  one  and  two  strip  computations  for  the 
analysis  of  the  flowfield  in  axisymmetric  nozzles  of  arbitrary  shape  with 
particular  emphasis  on  conical  nozzles.  The  two  strip  solutions  gave  good 
agreement  with  experimental  wall  and  centerline  pressure  measurements  for 
nozzles  with  R  >2.  Liddle  [72]  later  extended  the  method  to  include  the 
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cases  of  annular  and  asymmetric,  plane  nozzles.  One  strip  calculations 
were  performed  and  were  found  to  be  in  reasonable  agreement  with  wall  pres¬ 
sure  measurements  for  an  annular  configuration  and  with  the  numerical 
method  of  Katsanis  [73]  for  incompressible  flow  in  a  planar  nozzle.  One  of 
the  fundamental  difficulties  involved  with  using  the  method  of  integral 


relations  is  that  to  achieve  reasonable  accuracy  for  nozzles  which  have  a 
small  throat  radius  of  curvature  a  large  number  of  strips  is  required  with 
attendant  long  computing  times  and  program  complexity.  Saddle  point 
singularities  which  arise  in  the  integration  of  the  equations  must  also  be 
carefully  handled. 

Numerical  techniques  other  than  the  time  dependent  ones  discussed 
previously  have  also  been  utilized  for  throat  flowfield  analyses.  In  an 
early  investigation,  Emmons  [74]  developed  a  method  which  used  a  conformal 
transformation  and  a  relaxation  process  designed  for  hand  calculations  to 
solve  the  governing  equations.  In  a  later  publication  [75],  an  example 
was  calculated  for  the  flow  in  a  two-dimensional  nozzle  with  a  hyperbolic 
wall  contour.  Interestingly,  when  this  method  was  programmed  on  a  digital 
computer  convergence  problems  were  encountered  [15].  Twenty-five  years 
later  Prozan  and  Kooker  [76]  investigated  axisymmetric  nozzle  flowfields 
using  an  error  minimization  technique.  A  sophisticated  steepest  descent 
procedure  was  employed  to  minimize  the  residuals  in  the  steady  flow  con¬ 
servation  equations.  The  results  of  this  method  agree  quite  well  with  the 
flowfield  measurements  of  Cuffel ,  et  al.  [30]  for  a  nozzle  with  small 
throat  radius  of  curvature.  Typical  run  times  ranged  from  5  to  10  minutes 
on  an  IBM  7094.  Recently,  Brown,  et  al .  [77]  used  a  type-dependent 
relaxation  technique  to  study  both  conventional  and  annular  axisymmetric 
nozzle  flows  including  the  effects  of  rotational i ty.  In  this  method, 
which  is  quite  popular  for  external  transonic  flows,  the  axial  derivatives 
are  approximated  by  centered  differences  in  the  subsonic  (elliptic)  portion 
of  the  flowfield  and  by  backward  differences  in  the  supersonic  (hyperbolic) 


region.  The  results  of  example  calculations  for  hyperbolic,  conically  con¬ 
vergent,  and  annular  nozzles  indicate  that  the  rotational  effects  of 
entrance  flow  nonuniformities  can  be  significant  in  evaluating  nozzle 
performance.  The  example  problems  required  from  3h  to  5  minutes  and  a 
region  size  of  approximately  1 70K  on  an  IBM  370/158  computer.  Using  a 
somewhat  different  numerical  scheme,  Fanning  and  Mueller  [78]  have  studied 
planar  nozzle  flowfields.  The  problem  was  first  recast  as  a  Laplace  equa¬ 
tion  with  an  inhomogeneous  term  that  encompasses  the  compressibility  effects. 
A  successive  approximation  method  was  then  utilized  to  converge  to  a  final 
solution.  Flowfield  and  discharge  coefficient  results  agree  well  with 
experimental  measurements,  but  the  computing  times  are  extremely  long,  taking 
from  45  minutes  to  6  hours  of  IBM  370/155  time. 

B.  EXPERIMENTAL 

Not  only  are  theoretical  analyses  of  transonic  nozzle  flows  quite 
complicated,  but  experimental  measurement  of  them  is  also  difficult.  This  is 
because  disturbances  introduced  by  measurement  probes  are  propagated  along 
characteristics  and  are  reflected  from  the  duct  walls  back  into  the  region 
of  interest.  Thus,  most  experimental  studies  have  been  limited  to  wall 
static  pressure  measurements,  although  some  investigators  have  attempted  to 
use  pitot  probes.  Recently,  non-intrusive  optical  techniques  have  also  been 
employed. 

The  most  thorough  flowfield  experiments  have  been  performed  for  the 
two-dimensional  planar  geometry,  primarily  because  sidewall  static  taps  can 
be  used  to  obtain  both  transverse  and  streamwise  pressure  distributions. 
Jacobs  [79]  utilized  this  method  to  determine  isobars  in  six  asymmetric  plane 
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nozzles.  Later,  Flack  and  Thompson  [14,53,80]  conducted  an  extensive  series 
of  experiments  on  a  set  of  13  two-dimensional  nozzles  with  a  number  of  wall 
curvature  combinations.  For  each  configuration  the  sonic  line  location  was 
determined  using  sidewall  static  pressure  readings  and  also  laser  doppler 
velocimeter  (LDV)  measurements .  It  was  found  that  the  two  data  sets  and  the 
results  of  the  series  solution  developed  in  [37]  agreed  with  each  other 
reasonably  well  with  the  LDV  measurements  being  slightly  upstream  of  the 
static  pressure  measurements  and  the  theoretical  results  being  still  further 
upstream.  These  investigators  observed  that  small  changes  in  the  boundary 
geometry  caused  very  large  changes  in  the  sonic  line  location. 

Experiments  for  axisymmetric  Laval  nozzles  have  been  limited  mainly  to 
wall  static  pressure  readings  since  the  nozzle  geometry  does  not  provide 
any  means  for  making  flowfield  static  pressure  or  LDV  measurements .  One  of 
the  first  experimental  investigations  of  transonic  nozzle  flow  was  performed 
by  Stanton  [81].  The  flowfield  in  an  axisymmetric  nozzle  with  Rc  =  4.77 
was  studied  by  determining  static  pressure  distributions  along  the  contour 
for  a  number  of  supply  pressures.  Pitot  pressure  distributions  were  also 
measured  along  the  axis  and  along  an  axial  line  near  the  nozzle  wall,  but 
the  probe  certainly  disturbed  these  measurements,  particularly  those  made 
near  the  wall.  At  the  Jet  Propulsion  Laboratory  a  number  of  experimental 
investigations  [30,58,82-84]  of  axisymmetric  nozzle  flows  has  been  made  by 
Back  and  his  co-workers.  Static  wall  pressure  and  discharge  coefficient 
measurements  for  nozzles  with  a  dimensionless  throat  radius  of  curvature  as 
small  as  0.25  were  used  to  study  the  effects  of  contraction  section  shape, 
inlet  flow  direction,  wall  friction,  entrance  boundary  layer  thickness,  and 
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wall  cooling.  Only  for  a  nozzle  with  a  radius  of  curvature  of  =  0.625 
and  conical  inlet  and  exit  angles  of  45°  and  15°  ,  respectively ,  were  flow- 
field  measurements  made  [30],  These  were  obtained  by  means  of  a  0.889  mm 
diameter  tube  which  was  supported  and  held  taut  between  an  upstream  plenum 
chamber  and  a  downstream  vacuum  chamber.  The  tube  had  a  0.152  mm  diameter 
hole  drilled  radially  in  it,  and  the  pressure  distributions  along  the  axis 
and  along  five  radial  locations  were  found  by  appropriate  axial  traverses. 

At  the  off-axis  positions  disturbances  were  undoubtedly  introduced  into  the 
flowfield  and  the  velocity  vector  was  inclined  to  the  tube.  However,  to 
date  this  is  the  only  known  set  of  data  which  considers  the  flow  patterns 
in  the  transonic  region  of  an  axisymmetric  nozzle.  Liddle  and  Archer  [71] 
measured  the  centerline  and  wall  pressures  fo:"  four  axisymmetric  nozzles 
with  radius  of  curvature  values  ranging  from  0.625  to  4.0.  The  centerline 
distributions,  however,  are  reported  for  only  two  of  the  nozzles.  In 
another  series  of  experiments,  Masure,  et  al.  [63]  measured  the  discharge 
coefficients  for  five  nozzles  with  conical  inlet  and  exit  sections  and  R 

C 

values  from  0.25  to  0.8  . 

Data  for  annular  configurations  is  extremely  sparse.  Bresnahan  and 
Johns  [85]  measured  static  pressures  along  the  inner  wall  of  a  turbojet  plug 
nozzle.  In  their  experiments  the  outer  wall  was  straight  and  parallel  to 
the  axis  of  symmetry  and  the  inner  wall  had  a  large  radius  of  curvature, 

Rc  =  7.8  .  However,  there  were  only  three  pressure  taps  in  the  transonic 
region  for  which  the  corresponding  Mach  number  was  less  than  1.5,  since  the 
supersonic  region  was  of  primary  interest  in  this  study.  As  mentioned  pre¬ 
viously,  Brunell  [61]  conducted  experiments  on  an  annular  nozzle  with  a 


cy  1  indn'ca  1  centerbody,  but  only  operating  conditions  for  which  the  nozzle 
was  unchoked  were  considered.  Liddle  (72!  has  also  obtained  wall  pressure 
measurements  tor  an  annular,  supersonic  nozzle.  However,  the  geometry  of 
the  experimental  setup  is  not  specified  in  (721,  and,  in  any  event,  the  data 
are  very  close  to  the  results  of  one-dimensional  theory. 

Obviously,  there  is  a  need  to  expand  the  data  base  for  flow  in  the 
transonic  region  of  nozzles,  particularly  in  terms  of  flowfield  measurements 
for  both  conventional  and  annular  a xi symmetric  configurations.  As  mentioned 
in  the  preceding  chapter,  this  is  one  of  the  objectives  of  the  present 
investigation. 


III.  THEORETICAL  ANALYSIS 


The  objectives  outlined  in  the  introductory  chapter  for  the  theoretical 
investigation  are  to  provide  a  numerically  fast,  reliable,  and  accurate 
method  for  describing  the  throat  flowfields  in  annular,  supersonic  nozzles. 

In  addition,  it  is  desired  that  the  method  be  direct  and  that  it  be  appli¬ 
cable  to  a  widt  variety  of  configurations.  Of  the  few  previous  analyses 
which  have  considered  the  annular  geometry,  none  are  felt  to  have  all  of  the 
desired  characteristics .  Two  [23,25]  utilize  indirect  methods;  others 
[36,37,52,72]  may  be  applied  only  to  specialized  configurations,  and  the 
rest  [62,64,68,77]  are  generally  numerically  slow  and  expensive.  Only  the 
recent  numerical  techniques  of  Cline  [68]  and  Brown,  et  al.  [77]  have 
obtained  results  with  a  modest  expenditure  of  computing  time.  However,  both 
of  these  methods  implicitly  assume  that  the  main  flow  direction  in  the 
throat  is  essentially  parallel  to  the  axis  of  symmetry. 

Based  on  these  considerations,  the  method  which  has  been  chosen  is  a 
series  expansion  technique  similar  to  that  used  by  Moore  and  Hall  [36]  and 
Thompson  and  Flack  [37]  but  more  general  in  terms  of  the  geometries  which 
can  be  analyzed. 

A.  FORMULATION 

A  sketch  of  the  configuration  under  consideration  is  shown  in  Fig.  1 1 1 . 1 . 
The  R-Z  coordinate  system  is  the  standard  cylindrical  system  while  the  x-y 
coordinate  system  is  rotated  such  that  the  y-axis  lies  along  the  cross  section 
of  minimum  area  in  the  throat  with  the  x-axis  perpendicular  to  the  y-axis  and 
the  origin  located  on  the  axis  of  symmetry.  For  this  general  case,  the  main 


flow  /-direction  is  inclined  to  the  axis  of  symmetry  at  an  angle  ,  positive 
in  the  counterclockwise  direction,  and  d  is  the  throat  half-height  in  the 
R-Z  coordinates.  It  should  be  noted  that  for  an  inclined,  annular  nozzle 
the  minimum  area  cross  section  does  not  correspond  to  the  cross  section  of 
minimum  distance  between  the  contours.  The  influence  of  the  radial  coordi¬ 
nate,  R,  in  computing  the  annular  area  results  in  the  minimum  area  section 
being  nearer  the  axis  of  symmetry  than  the  minimum  distance  section. 

It  is  assumed  that  the  problem  to  be  analyzed  is  the  steady,  inviscid, 
i rrotational ,  adiabatic  flow  of  a  perfect  gas  in  an  annular  nozzle  operat¬ 
ing  in  the  choked,  continuously  accelerating,  Meyer  regime.  These  are  the 
usual  assumptions  and  operating  conditions  of  interest.  The  inviscid 
assumption,  in  particular,  is  justified  since  the  Reynolds  number  based  on 
the  throat  height  is  large  for  most  practical  applications  and  the  stream- 
wise  pressure  gradient  is  favorable.  Therefore,  the  boundary  layers  in  the 
throat  region  are  very  thin  and  viscous  effects  may  be  safely  neglected. 

Under  these  assumptions,  the  governing  equations  may  be  taken  as  the 
i rrotational i ty  condition  and  the  "gas  dynamic  equation"  [86], 

UR  -  vz  =  0  (III-l) 

(u2-a2)uz  +  (v2-a2)vR  +  2UVUR  -  ^  =  0  (HI-2) 

a  =  speed  of  sound  , 

where  the  gas  dynamic  equation  is  obtained  from  the  continuity  equation, 
Euler's  momentum  equation,  and  the  definition  of  the  speed  of  sound.  In 
these  equations,  subscripts  are  used  to  denote  partial  differentiation  with 
respect  to  Z  and  R.  These  two  governing  equations  can  be  combined  to  form 


u  =  1  +  u 


(III-10) 
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v  =  v  ,  (III-ll) 

where  u  and  v  are  assumed  to  be  small  compared  to  unity,  jui,  ivi^l.  This 
results  in 

u  -  v  =  0  (IU-12) 
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The  boundary  conditions  in  this  inviscid  analysis  are  that  the  bounding 
walls  must  be  streamlines.  Using  y  =  g(x)  and  y  =  h(x)  as  the  functional 
forms  of  the  equations  for  the  inner  and  outer  wall  contours,  respectively, 
the  boundary  conditions  can  be  stated  as 

v(x,g(x) )  =  [l+u(x,g(x) )] g*  (x)  (IH-14) 

v (x ,h(x ) )  =  [l+u(x,h(x))]h‘ (x)  ,  (III-15) 

where  primed  quantities  represent  differentiation  with  respect  to  x, 

(  )'  =  ^  .  Expanding  the  equations  for  the  boundaries  in  Maclaurin  series 
about  the  throat. 


2  3 

y  =  g(x)  =  y.  +  g'(0)x  +  g"(0)  |r  +  g"'  (o)  J r  +  ... 

2  3 

y  =  h ( x )  =  yQ  +  h'(0)x  +  h"(0)  | r  +  h,H  (0)  |r  +  . . . 
the  boundary  conditions  can  be  rewritten  as 

v(x,g(x) )  =  [l+u(x,g(x))]  j^g '  (0)+g" (0)x+g'M  (0)  |y  + 
v(x,h(x))  =  [l+u(x,h(x) )]  |^h '  (0)+h"(0)x+h'"  (0)  jr  + 


( 1 1 1-16) 
( 1 1 1-17 ) 


( II 1-18) 
(III-19) 
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The  solution  technique  involves  expanding  the  perturbation  velocity 
components  in  appropriate  series  and  substituting  them  into  the  governing 
equations  ana  boundary  conditions.  However,  before  this  can  be  done  an 
expansion  parameter,  t ,  must  be  chosen  and  various  order  of  magnitude  esti¬ 
mates  made.  Both  Kliegel  and  Levine  [46]  and  Thompson  and  Flack  [37]  found 
that  improved  convergence  properties  of  their  series  solutions  were  attained 
for  nozzles  with  a  small  wall  radius  of  curvature  by  using  =  (Rc+1)  1 
rather  than  the  Hall  [35]  parameter  c  =  R.  '.  Based  on  their  experience, 
the  expansion  parameter  which  has  been  used  in  this  investigation  is 

r  =  (R  +n)_1  ,  ( II 1-20) 
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where  Rc  is  an  average  dimensionless  radius  of  curvature  for  the  two  bounding 
walls  and  n  is  arbitrary.  For  r,  >  1,  e  is  therefore  less  than  unity  regard¬ 
less  of  how  small  R^  may  be.  Defining  Rc  in  terms  of  the  second  derivatives 
of  the  equations  for  the  wall  contours, 

*.  =  h"  (0)  -  g"  (0)  •  (  m-21 ) 

t:  becomes 


£  = 


-Jim  -..gloj. 

+  n[h"(0)  -  g"(0)] 


( II 1-22) 


This  is  identical  to  the  parameter  used  in  [37]  except  that  n  may  have  any 
value  here  and  is  not  restricted  to  the  values  zero  or  one. 

From  one-dimensional  considerations,  it  is  well  known  [13]  that  the 
acceleration  at  the  throat  of  a  supersonic  nozzle  is  given  by 

I  1/2 
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X  =  Q 

(y+1)Rc 

( 1 11-23) 
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where  j  =  0  for  plane  two-dimensional  nozzles 
=  1  for  axisymmetric  nozzles 

Lis  ng  La.  (III-10),  the  definition  of  transonic  perturbation  velocity  com¬ 
ponent  u ,  and  Eq.  (Ill- 20 ) ,  t  .e  definition  of  the  expansion  parameter  e, 
Eq.  (III-23)  yields  the  following  order  of  magnitude  estimate  for  the 
transonic  throat  region. 


Also  from  the  dominant  terms  in  boundary  conditions  (III-18)  and  (III-19),  it 
is  known  that 

vk  |  =  9"(0)  +  g1"  (0)x  +  ...  (II 1-25 ) 

vjy=h  =  h"  (0)  +  h"'  (0)x  +  ...  .  ( 1 1 1-26) 

It  is  shown  in  Appendix  A  that  the  terms  g"(0)  and  h"(0)  dominaie  the  right 
hand  sides  of  the  two  preceding  equations  for  nozzles  with  circular  arc  wall 
contours.  Combining  this  information  with  the  definition  of  e  in  Ea.  ( 1 1 1-22 ) 
allows  the  estimate 

v  =  0(e)  ( 1 1 1-27 ) 

X 

for  the  transonic  region. 

Requiring  that  both  terms  in  the  irrotational ity  condition  (III-12)  be 
of  the  same  order  of  magnitude  and  using  this  last  estimate  for  vx  results 
in  the  conclusion  that  Ciy  must  also  be  0(e)  in  the  throat.  Since  all  lengths 
have  been  non-dimensionalized  with  respect  to  the  throat  half-height,  d.  Ay 
is  o(l)  across  the  throat  and  therefore  u  must  be  0(e)  in  the  transonic  region 
of  interest.  Utilizing  the  estimates  expressed  in  Eqs.  (II 1-24)  and  (II 1-27 ) , 
it  is  concluded  that  v/u  =  o(eI|/2)  or  v  =  o(e3/,2j.  Taking  this  estimation  for 
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the  magnitude  of  v  together  with  Eq.  (III-27),  or  alternately  the  irrotation- 
ality  governing  equation  (III-12),  results  in  the  requirement  that  x  be  of 
-  [i*/2]  for  consistency.  Thus,  the  solution  developed  here  is  restricted  to 
the  narrow  transonic  region  about  the  throat  plane,  x  =  0. 

To  complete  the  order  of  magnitude  analysis  the  gas  dynamic  equation 
(1 11-13)  must  be  investigated.  Dividing  the  numerator  and  denominator  of 
the  last  term  of  this  equation  by  cosB,  it  may  be  restated  as. 


cosB  t  0  .  ( 1 11-28) 

Disregarding  their  numerical  coefficients,  the  dominant  factors  in  each  of 
the  first  three  terms  are  uu  ,  vu  ,  and  v  which  are  c? (e3 and 

x  y  y 

(r.3//2j,  respectively.  Likewise,  the  dominant  factors  in  the  last  term  are, 

v  +  tang 
y  +  x  tang  ' 

Assuming  for  the  moment  that  tang  is  at  most  0(1),  i.e.,  disallowing  con¬ 
figurations  with  radial  or  nearly  radial  main  flow  directions  for  which 
tan8  +  oo*  ,  this  term  can  further  be  reduced  to 

v  +  tanB 

y 

since  x  has  been  restricted  to  the  throat  region,  x  =  ofc1^2).  Since 
y  =  0(1)  or  larger,  depending  on  the  dimensionless  distance  from  the  axis  of 

fThis  restriction  is  consistent  with  the  previous  one  that  cosB  f  0  in 
Eq.  ( 1 1 1-28) . 
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symmetry  to  the  annular  nozzle  throat,  and  using  the  previous  estimate  for 


the  magnitude  of  v,  it  i 


or  smaller.  Therefore, 
be  of  the  same  order,  it 
tanB  =  ^  [e3/2 ) 


known  that, 

n  order  that  the  dominant 
is  required  that 


( 1 1 1-29) 

terms  in  Eq.  ( 1 1 1-28)  all 

( 1 1 1-30) 


or  smaller.  This  means  that  for  annular  nozzles  with  throats  that  are  near 
the  axis  of  symmetry  in  a  dimensionless  sense,  y  =  o(l),  only  small  angles 
of  inclination,  6,  may  be  considered.  However,  this  would  seem  to  be  the 
only  physically  realistic  case  anyway,  since  one  would  not  expect  to 
encounter  an  annular  configuration  for  which  the  throat  is  very  close  to 
the  symmetry  axis  while  the  main  flow  x-direction  is  highly  inclined  to  it. 
For  nozzles  whose  throats  are  large  distances  from  the  axis,  y  »  1,  the 
restriction  to  small  angles  of  inclination  may  be  relaxed  as  long  as  the 
estimate  of  Eq.  (II 1-30 )  is  satisfied.  In  the  limit  as  y  -*■  °°,  the  last  term 
in  the  governing  gas  dynamic  equation  (II 1-28)  vanishes  as  it  takes  the 
correct  planar  form. 

The  expansion  parameter,  e,  as  defined  in  Eq.  (II 1-22 ) ,  is  related  to 
the  average  wall  curvature  of  the  bounding  walls  at  the  throat  and  enters 
the  formulation  through  the  second  derivative  terms  in  the  Maclaurin  series 
expansions  for  the  boundary  conditions,  Eqs.  ( 1 1 1-18)  and  (III-19).  However, 
because  of  the  complication  that  the  definition  of  e  involves  both  g"(0)  and 
h"(0)  while  each  boundary  condition  contains  only  one  or  the  other,  a  method 
must  be  devised  such  that  when  expansions  are  assumed  for  the  dependent 
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variables,  appropriate  series  occur  in  the  expression  of  the  boundary  condi¬ 
tions.  This  is  done  in  the  following  manner.  Solving  Eq.  (Ill- 20 )  for  R 


and  using  the  identity  g"(0)  =  g“ ( 0 )  Rc  Rc  1  along  with  Eq.  ( 1 1 1 - 2 1 )  for  R 
and  Eq.  (III-31)  for  R^  1  yields 


2gn(0) 

h“ (0)  -  g" (0) 


2  2  3 

ie+ne  +n  e  + 


■  j,  ( I  i  1-32) 


for  ) ne |  <1  . 


Simi 1 arly. 


h"(0)  = 


2h" (0)  £ 

PIT  g"(0) 


2h" (0) 

nc  "  h" ( 0 )  -  g" (0) 


,  2  2  3 

(c+ne  +n  c  + 


( 1 1 1-33) 


for  | ne |  <  1  . 

The  requirement  that  |ne|  <  1  in  the  expansions  above  is  always  met  since 
h"(0)  -  g"(0)  >  0  for  an  area  minimum  at  the  throat  and  only  non-negative 
values  of  n  are  considered  (see  Eq.  (III-22)).  In  this  way  power  series  in 
€  are  obtained  in  the  boundary  condition  evaluations,  although  it  is  clear 
that  the  second  derivatives  are  exactly  matched  only  for  n  =  0.  The  coeffi¬ 
cients  multiplying  the  parentheses  in  Eqs.  ( 1 1 1-32 )  and  (II 1-33 )  are  0(1), 
and  as  shown  in  Appendix  A,  the  first  few  derivatives  in  the  Maclaurin 
series  expansions  for  circular  arc  wall  contours  are  of  the  following  orders 
of  magnitude, 

g'(0),  h‘ (0)  =  o(e3/2)  g"(0),  h"(0)  =  o(e) 

g'"  (0),  h'"  (0)  =  o[t/2)  giv  (0),  hiv  (0)  =  o(e3) 


(III- 34 ) 


Note  particularly  the  transformation  from  coordinate  x  to  the  0(1),  stretched 
coordinate  z.  It  is  also  seen  that  the  estimates  y  =  0(1)  and  tang  =  tf(e  3^2 
have  been  used  in  the  definitions  of  y  and  ,  although  as  previously  dis¬ 
cussed,  nozzles  with  larger  inclinations  may  be  analyzed  if  the  requirement 
expressed  in  Eq.  (III-30)  is  satisfied. 

With  these  definitions  and  order  of  magnitude  estimates  the  next  step  is 
to  expand  the  perturbation  velocity  components  in  appropriate  series. 

Assuming  expansions  in  half-powers  of  e  starting  with  u  =  0(e)  and 
v  =  4c5/2),  it  has  been  found  that  the  following  series  are  sufficient, 

u(z,y)  =  Uj(z,y)c  +  u2(z,y)e2  +  u3(z,y)e3  +  ...  ( 1 1 1-36 ) 

-  .  -11/2  r  1 

?(z,y)  =  e  VjU.yJe  +  v2(z,y)e2  +  v3(z,y)e3  +  ...  ,  ( 1 1 1-37 ) 

i.e.,  terms  involving  odd  half-powers  of  e  in  the  u  expansion  and  even  half¬ 
powers  of  e  in  the  v  expansion  vanish. 
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Substituting  definitions  ( II 1-35)  and  expansions  (III- 36 )  and  ( 1 1 1 - 37 ) 
into  the  governing  irrotational ity  equation  (III-12)  and  gas  dynamic  equa¬ 
tion  (II 1-28)  and  gathering  coefficients  of  like  powers  of  •  results  in  the 
following  two  sets  of  equations, 

3u  3v 

-£■  -  -gj-  =  0  (n  =  1,2,3,  ...)  ( 1 1 1-38) 

and 


3u  3v  6  +  v 

'  2“i  ir  +  17  +  -7~ 


(n  =  1) 


(111-39) 


3u  3u  3v  v 

-  2u,  —l1-  -  2u  — -  +  — T-— ■  +  —  =f  fu  ,v,  ,  ...,u  ,v  )  (n  =  2,3,  . .  .  ) , 

1  dZ  n  3Z  3y  y  n  (  1  ’  1  n-l  n-1  j 

(1 11-40) 

where  the  functions,  f  ,  on  the  right  hand  side  of  Eq.  (Ill -40 )  are  always 

n 

known  from  the  lower  order  solutions.  The  first  two  are  given  by 
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In  the  last  term  of  gas  dynamic  equation  (II 1-28) ,  the  binomial  expansion  of 
the  factor  (y+x  tang)’1  has  been  used  which  assumes  that  |x  tanB/y|  <  1. 
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This  requirement  is  certainly  met  if  the  one  stated  in  Eq.  (Ill- 30 )  is.  Of 
interest  is  the  fact  that  the  irrotational ity  condition  retains  its  simple, 
linear  form  through  all  orders  and  that  the  only  nonlinear  equation  among 
the  set  of  gas  dynamic  governing  equations  is  the  first  order  one.  This 
behavior  is  typical  of  expansion  solution  techniques. 

The  boundary  conditions  which  apply  to  each  solution  order  must  also  be 
obtained  in  an  appropriate  form.  Only  the  details  for  the  boundary  condi¬ 
tion  at  the  inner  wall  will  be  given  since  the  procedure  for  the  outer  wall 
boundary  condition  is  identical.  Transforming  from  coordinate  x  to  coordi¬ 
nate  z,  Eq.  (III-14)  for  the  inner  wall  condition  becomes, 

v(z,g(z))  =  [1  +  u(z,g(z))]  |jf  ( 1 11-43) 

where  the  Maclaurin  series  expansion  for  the  inner  contour,  Eq.  ( 1 1 1 - 1 6 ) ,  is 
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(II 1-44) 


when  the  definitions  and  expansions  in  Eqs.  (II 1-32) ,  (II 1-33 ) ,  and  (II 1-35) 
are  introduced.  The  derivative  dg/dz  needed  in  Eq.  (III-43)  is. 
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Expansion  of  the  perturbation  velocity  components  at  the  inner  wall  in  Taylor 
series  about  y  =  y.  yields. 
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and 
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The  procedure  then  is  to  substitute  these  expansions  into  boundary  condition 
( 1 11-43)  using  the  series  for  u  and  v,  Eqs.  (Ill- 36 )  and  (III- 37 ) ,  and  the 
wall  contour  information,  Eqs.  (Ill- 44 )  and  (III-45).  Gathering  and  equat¬ 
ing  coefficients  of  powers  of  t.  leads  to  the  following  set  of  inner  wall 
boundary  conditions, 
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The  boundary  conditions  along  the  outer  contour  are  exactly  analogous, 

V1  (Z’yo)  =  hl  +  h2Z  (III-51) 

V2(Z**o)  =  h2nZ  +  (hl+h2Z)Ul  (Z’-0  (  I  I  1-52  ) 
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In  performing  these  Taylor  series  expansions  about  the  inner  and  outer  throat 
wall  locations,  it  is  clear  that  the  boundary  conditions  are  satisfied  exactly 
only  at  these  points. 
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The  formulations  for  the  various  orders  of  solution  in  the  expansion 
technique  are  now  complete.  For  the  first  order,  n  =  1,  irrotational i ty 
condition  (111-38),  gas  dynamic  equation  (111-39),  and  boundary  conditions 
(111-48)  and  (III-51)  comprise  the  formulation.  Likewise,  Eqs.  (Ill- 38 ) , 
(III- 40 ) ,  (III- 49 ) ,  and  (III-52)  are  the  ones  of  importance  for  the  second 
order,  etc.  .  Once  the  solutions  for  the  various  orders  in  the  expansion 
technique  have  been  obtained,  the  flowfields  in  the  throat  region  of 
annular,  supersonic  nozzles  may  be  characterized. 


B.  SOLUTION 

The  solution  technique  is  similar  to  that  used  originally  by  Hall  [35] 
and  proceeds  by  first  considering  the  formulation  for  the  first  order  prob¬ 
lem  and  then  moving  to  the  second  and  higher  order  problems.  As  just  men¬ 
tioned,  the  governing  equations  and  boundary  conditions  for  the  first  order 
are 
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vi (z’yi  ]  =  +  92z  ( 1 11-48) 

vi  (*.*.)  =  h!  +  V  •  <ni-51> 

Inspection  of  boundary  conditions  (II 1-48)  and  ( 1 1 1-51 )  suggests  the 
following  solution  form  for  v  , 


Vj (z,y)  =  A^(y)  +  A ^ (y )z  , 


(III- 54 ) 


where 
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and  the  prime  is  now  used  to  denote  differentiation  with  respect  to  y. 
Substitution  of  the  v  solution  form  into  i rrotationa 1 i ty  equation  (III- 38 ) 
and  integration  with  respect  to  y  results  in 

u,(z,y)  =  Aj(y)  +  BQ  +  BlZ  ,  ( 1 11-57) 

where  BQ  and  8  are  constants  and  the  form  of  the  integrating  function  is 
suggested  by  Eq.  (Ill- 54 )  for  v  .  The  assumption  of  terms  involving  higher 
powers  of  z  in  either  Eqs.  (Ill- 54 )  or  (III- 57 )  is  unnecessary  since  these 
terms  would  vanish  by  boundary  conditions  (III -48 )  and  (III-51)  anyway. 

The  unknown  A  functions  and  B  constants  must  still  be  determined.  This 
is  accomplished  by  first  substituting  the  expressions  for  Uj  and  v(  ,  into 
the  gas  dynamic  governing  equation,  (III-39),  and  gathering  coefficients  of 
similar  powers  of  z.  Since  the  resulting  equation  is  to  be  satisfied  for  all 
z,  these  coefficients  must  vanish  which  leads  to  the  following  two  equations, 

(1 11-58) 
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In  this  manner  the  problem  has  been  reduced  from  the  solution  of  the  original 
set  of  partial  differential  equations  and  boundary  conditions  to  the  solution 
of  the  two  preceding  ordinary  differential  equations  with  conditions  (III- 55 ) 
and  (II 1-56 ) .  Equation  (II 1-58)  can  be  integrated  once  to  give 
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(II 1-60 ) 
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where  the  constants  B2  and  B  are  determined  from  the  conditions  in 
Eq.  ( 1 1 1-55)-  as 
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The  positive  root  is  chosen  in  the  expression  for  B  since  Su/3z,  which  is 
identically  B( c  for  the  first  order  solution,  must  be  positive  for  an 
accelerating  nozzle  flow.  Investigation  of  the  radicand  reveals  that  it 
also  is  always  positive  for  the  nozzle  configurations  considered  here,  h2  >  0 
and  g2  <  0.  Integration  of  Eq.  (Ill- 60 )  yields 

A,  (y)  =  \  B^y2  +  B2lny  ,  (III-63) 

where  an  integration  constant  is  unnecessary  since  it  can  be  considered  as 
being  included  with  the  still  undetermined  constant  Bq  in  Eq.  ( 1 1 1-57 )  for  ut  . 

Substituting  this  relation  for  (y)  into  Eq.  (Ill- 59 )  and  integrating 
once  yields 
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The  conditions  listed  in  (III- 56)  may  then  be  used  to  determine  B 
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where  B4  and  B5  have  been  defined  as. 
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The  first  order  solution  is  complete.  The  B  constants  are  evaluated  in 
the  following  order:  B2  ,  Eq.  (II 1-61 ) ;  B(  ,  Eq.  (Ill -62 ) ;  B4  ,  Eq.  (III-67); 
B&  ,  Eq.  (Ill- 68 ) ;  B3  ,  Eq.  (Ill- 65 ) ;  and  BQ  ,  Eq.  (II 1-66 ) .  With  these  con¬ 
stants  determined,  the  functions  (y),  A t(y),  and  A^(y)  are  then  known  from 
Eqs.  (Ill- 60 ) ,  (II 1-63 ) ,  and  (II 1-64 ) ,  and  ut  and  vj  may  be  evaluated  from 
Eqs.  ( 1 1 1- 57 )  and  (III-54),  respectively. 

The  higher  order  solutions  are  obtained  in  a  very  similar  manner  start¬ 
ing  from  the  formulations  developed  in  the  preceding  section.  Solutions 
through  the  first  three  orders  have  been  found,  and  as  might  be  guessed  from 
the  form  of  the  gas  dynamic  governing  equations,  Eqs.  (  1 1 1 - 39 ) - (  1 1 1-42) ,  the 
effort  required  to  obtain  them  increases  dramatically  as  one  moves  to  the 
higher  order  problems.  For  the  second  order,  v2  and  u2  are  of  the  form, 

v2(z,y)  =  C^y)  +  C;(y)z  +  C’(y)z2  (III-69) 


u2(z,y)  =  jC,  (y)+D0j 
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requiring  the  solution  of  three  ordinary  differential  equations,  while  the 
third  order  problem  requires  solving  four  ordinary  differential  equations 
with  v3  and  u3  given  by, 
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v3(z,y)  =  E‘(y)  +  E;(y)z  +  E^(y)z2  +  (y)z3 


( 1 1 1- 71 ) 


u  (z,y)  =  |E  (y)+F0i  +  |2E  (y)+F,'z  +  3E3(y)+F2  z2  +  F3z3  .  (111-72) 

Lengthy  algebraic  manipulations  are  involved,  but  the  form  of  the  differential 
equations  to  be  solved  is  very  similar  to  Eqs.  (Ill- 58 )  and  (III- 59 ) .  As  in 
the  integration  of  the  equations  for  the  first  order  problem,  a  large 
number  of  intermediate  constants  are  defined  and  used,  but  because  of  space 
limitations  the  details  of  the  second  and  third  order  solutions  will  not  be 
presented  here.  However,  a  complete  solution  summary  is  included  in  Appendix  B. 
With  the  expressions  for  the  fu^vl,  (u2,v2j,  and  fu3,v3j  perturbation 


velocity  components  determined,  other  quantities  of  interest  may  also  be 
found.  These  include:  the  velocities  components  u  and  v  in  the  x-y  coordi¬ 
nate  system;  M*,  the  ratio  of  the  local  speed  to  the  critical  speed  of  sound; 
6,  the  angle  of  inclination  of  the  velocity  vector  from  the  x-axis;  the  Mach 
number,  M;  and  the  local  static  to  stagnation  pressure  ratio,  P/P0 •  The 
series  expansions  for  these  quantities  in  terms  of  the  perturbation  velocity 
components  are  given  below, 
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The  reason  that  these  expansions  should  be  employed  is  so  that  the  only  terms 

included  are  those  for  which  the  coefficients  of  the  various  powers  of  e  are 

completely  known.  For  example,  if  u  and  v  are  first  calculated  from 

Eqs.  ( 1 1 1 - 73 )  and  ( 1 1 1- 74 )  and  then  M*  is  determined  from  the  definition 

f  2  2]1/2 

M*  =  [u  +v  J  ,  information  from  just  one  of  the  many  sixth  order  terms. 


VjC6 ,  will  be  included. 

Another  quantity  of  importance  is  the  discharge  or  flow  coefficient,  CD , 
defined  as  the  ratio  of  the  actual  nozzle  mass  flow  to  that  obtained  from  the 
ideal  one-dimensional  approximation  of  uniform,  sonic  flow  at  the  throat,  i.e.. 
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where  in  the  present  notation  u  has  been  non-dimensional ized  with  respect  to 
the  critical  speed  of  sound,  a*.  Substituting  the  series  expansions  for  p/p* 
and  u  and  the  expression  for  dA/A*  into  the  integrand,  the  relation  for  CD 
becomes 
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after  integration,  where  the  constants  C  ,  CD2 ,  and  CDJ  are  evaluated  as 
shown  in  Appendix  B.  The  reduction  in  mass  flow  due  to  the  two-dimensional, 
nonuniform  nature  of  the  flowfield  is  generally  small  and  the  value  of  the 
discharge  coefficient  is  therefore  very  close  to  unity  for  the  configurations 
which  may  be  analyzed  with  the  expansion  solution  developed  here. 

Thus,  with  the  fu^v^,  |^2*v2j»  and  ^u3’v3j  transonic  perturbation 
velocity  components  and  the  C  ,  CD2 ,  and  CD3  discharge  coefficient  constants 
determined,  all  of  the  flow  variables  of  interest,  including  those  just  pre¬ 
sented,  are  known  to  the  third  order  in  the  present  series  approximations. 

C.  SOLUTION  CHECKS 

Because  the  development  of  the  solutions  described  in  the  preceding 
section  requires  lengthy  algebraic  manipulations,  particularly  for  the  third 
order  solution,  they  must  be  thoroughly  checked  to  ensure  their  validity. 

To  accomplish  this  task  two  groups  of  tests  have  been  carried  out  as  shown 
in  Table  III.l. 

The  first  set  of  tests  involves  reducing  the  solution  obtained  here  for 
the  general  annular  configuration  to  previous  solutions  for  simpler 
geometries.  The  results  for  the  conventional  axisymmetric  configuration  are 
found  by  passing  to  the  limit  y.  -*  0  for  a  nozzle  with  a  straight  inner 
boundary,  i.e.,  as  the  inner  wall  approaches  the  axis  of  symmetry.  The 
limiting  value,  y(  =  0,  cannot  be  used  in  the  numerical  implementation  of 
the  annular  solution  since  this  leads  to  division  by  zero.  However,  an 
arbitrarily  small  value,  e.g.,  y;  =  10”10,  may  be  employed  to  approximate  the 
axisymmetric  configuration.  In  a  similar  manner,  the  planar  symmetric 
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geometry  may  be  investigated  by  considering  the  limit  y  *  since  the 
transverse  curvature  effect  becomes  negligible  in  that  limit.  In  this  sit¬ 
uation  there  are  two  possible  wall  configurations  of  interest.  Both  the 
case  in  which  the  bounding  walls  have  equal  and  opposite  curvatures  and 
that  in  which  one  of  the  walls  is  straight  approach  the  symmetric  plane 
limit  as  y.  'c;  in  the  latter  case,  the  straight  wall  may  be  considered  as 
the  centerline  of  a  symmetric  nozzle.  Using  n  =  0  in  the  expansion  param¬ 
eter  definition,  (Eq.  1 1 1-22) .  the  results  of  the  present  solution  reduce  to 
those  of  Hall  [35],  with  corrections  by  Kliegel  and  Levine  [46],  in  both 
the  ax i symmetric  and  plane  symmetric  limits.  Likewise,  for  n  =  2  the  solu¬ 
tion  developed  here  is  identical  to  previous,  unpublished  solutions  [87]  by 
the  present  author  for  both  the  conventional  axisymmetric  and  plane  configur¬ 
ations,  where  e  =  jRc+lj_1  was  employed  as  the  expansion  parameter.  It 
should  be  emphasized  that  neither  of  these  e  =  ^Rc+lj_1  solutions  for  the 
axisymmetric  case  coincide  with  that  of  Kliegel  and  Levine  [46]  because,  as 
pointed  out  in  Chapter  II,  their  result  is  not  a  solution  since  it  doesn’t 
satisfy  the  equations  of  motion.  Rather,  it  is  simply  a  transformation  from 
the  Hall  series  in  e  =  Rc_1  to  a  series  in  e  =  j^+lj'1  such  that  the  two 
are  equivalent  as  Rc  ♦  ».  Except  for  errors  in  the  third  order  discharge 
coefficient  constants,  the  analysis  of  Thompson  and  Flack  [37]  for  plane 
nozzles  can  also  be  reduced  to  the  e  =  ^Rc+lj_1  so^utions  mentioned  above 
for  the  symmetric,  planar  case.  Another  configuration  which  provides  a  limit 
check  on  the  solution  for  the  general  annular  nozzle  is  the  specialized  case 
for  which  either  the  inner  or  outer  boundary  is  straight  and  parallel  to  the 
axis  of  symmetry.  Results  of  the  two  term  solutions  developed  in  [87]  using 
e  =  [Rc+l|_1  agree  with  those  of  the  present  solution  for  n  =  2. 
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It  is  to  be  noted  that  when  one  of  the  bounding  walls  is  straight,  e.g., 
g"(0)  =  0,  the  expression  for  c,  Eq.  (Ill - 22 ) ,  simplifies  to 


b"  (0) 

2  +  nh"(0)  ’ 


( 1 1 1-81 ) 


which  is  equivalent  to 

t  =  [2Rc+n]_1  ,  (II 1-82) 

since  h"(0)  =  R  -I  if  h'(0)  =  0.  Because  of  the  presence  of  the  factor  2  in 

C 

Eq.  (Ill- 82 )  it  might  be  expected  that  for  cases  with  one  straight  boundary, 
different  results  would  be  obtained  with  the  present  analysis  than  with  pre¬ 
vious  analyses  which  use  either  e  =  Rc_1  or  e  =  |rc+iJ~'  as  the  expansion 
parameter.  However,  if  the  formulation  of  the  present  investigation  is  care¬ 
fully  studied,  it  is  found  that  n  enters  the  problem  only  as  the  product  ge 
through  the  boundary  conditions,  Eqs.  (III-32)  and  (III-33).  As  long  as  this 
product  is  the  same  for  two  solutions,  the  results  obtained  with  them  will  be 
identical  regardless  of  the  numerical  coefficient  of  Rc  in  the  expansion 
parameter.  Therefore,  no  generality  is  lost  by  employing  the  definitions 
given  in  Eqs.  ( 1 1 1-20)- ( II 1-22)  for  e.  This  explains  why  all  of  the  n  =  0 
solutions  produce  identical  results  independent  of  the  e  definition  and  also 
why  the  n  =  2  case  in  the  present  technique  corresponds  with  the  previous 
e  =  |Rc+lj_1  analyses. 

The  reduction  tests  just  discussed  were  generally  carried  out  numerically 
by  comparing  the  discharge  coefficient  constants  CD1  ,  CD2 ,  and  CD3  and  the 
transonic  perturbation  velocity  components  [uj.Vjj,  ^u2,v2j,  and  (u3»v3] 
obtained  with  the  present  solution  to  those  obtained  from  earlier  solutions 
for  the  specialized  configurations.  As  demonstrated  by  the  expansions  given 
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in  Eqs.  (  1  1 1-73)- (  1  1  1-78)  and  (III- 80 ) ,  these  parameters  completely  specify 
the  quantities  of  interest.  The  perturbation  velocity  comparisons  were 
made  by  evaluating  the  components  at  a  number  of  y  locations  on  planes  of 
constant  x  through  the  throat  region.  Only  for  the  limiting  case  of  the 
conventional  axisymmetric  nozzle  is  it  convenient  to  simplify  the  general 
solution  analytically.  Since  this  is  a  special  case  of  great  practical 
importance,  the  resulting  solution  is  summarized  in  Appendix  C,  where  the 
expansion  parameter  has  been  taken  as  e  =  ^c+'ij_1-  For  n  f  0  this  is  a  new 
and  very  useful  result. 

The  second  group  of  tests  involves  numerically  back  substituting  the 
solutions  for  the  various  orders  .into  the  corresponding  governing  equations, 

( 1 1 1 - 38) - ( 1 1 1 -42 ) ,  and  boundary  conditions,  ( 1 1 1-48)- ( 1 1 1-53)  and  evaluating 
the  residuals.  The  partial  derivatives  in  these  equations  are  approximated 
by  using  second  order,  central,  finite  differences  for  all  z-derivatives  and 
for  y-deri vati ves  at  interior  points  and  either  forward  or  backward,  second 
order  differences  for  the  y-deri vatives  at  boundary  points.  This  scheme  pro¬ 
vides  a  powerful  means  of  checking  the  solution  given  its  formulation  and  can 
be  used  to  verify  it  in  this  sense. 

The  satisfaction  of  the  governing  equations,  (II 1-38)- ( III-42) ,  as 
measured  by  the  size  of  the  residuals,  was  tested  at  a  number  of  y-locations 
ranging  from  the  inner  to  the  outer  boundary  on  planes  of  constant  x  through 
the  throat  region.  The  boundary  conditions  were  also  checked  at  various  x 
stations  through  this  region.  In  addition,  the  tests  were  carried  out  for  a 
large  number  of  geometrical  configurations  and  specific  heat  ratios  for  the 
gas  so  that  all  of  the  significant  parameters  were  varied.  In  all  cases  it 
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was  found  that  the  residuals  for  the  equations  and  boundary  conditions  were 

f  -8  I  -4 

0(10  j  or  smaller  when  a  grid  spacing  of  J.x  =  Ly  =  10  was  used  for  the 
di fferenci ng. 

The  results  of  the  two  sets  of  tests  described  above  provide  strong 
evidence  of  the  correctness  of  the  solution  developed  in  the  preceding  sec¬ 
tion.  It  is  felt,  therefore,  that  the  results  of  this  analysis  can  be  used 
with  confidence  to  predict  the  flowfields  in  the  throat  regions  of  annular, 
supersonic  nozzles. 

D.  PARAMETRIC  STUDY  OF  SOLUTION  BEHAVIOR 

So  far  no  mention  has  been  made  either  of  the  convergence  properties  of 
the  series  solution  which  has  been  developed  or  its  range  of  applicability 
to  various  geometrical  configurations,  etc.  .  The  purpose  of  this  section 
is  to  present  the  results  of  a  comprehensive,  numerical  study  whose  objec¬ 
tive  was  to  parametrically  investigate  these  questions. 

As  mentioned  in  the  first  section  of  this  chapter,  the  order  of  magni¬ 
tude  estimate,  y  =  <9(1 )  or  larger,  was  made,  and  the  origin  of  the  x-y 
coordinate  system  was  taken  to  lie  on  the  axis  of  symmetry  so  that  annular 
configurations  whose  throats  are  near  the  axis  can  be  analyzed.  These 
geometries  are  of  prime  importance  for  the  supersonic-supersonic  ejector 
application,  and  use  of  this  estimate  also  means  that  the  solution  may  be 
applied  to  conventional  axisymmetric  nozzles.  However,  for  annular  axi- 
symmetric,  or  in  the  limiting  case,  planar  configurations  for  which  the 
dimensionless  distance  from  the  symmetry  axis  to  the  throat  region  is  large, 
y  »  1,  the  evaluation  of  the  transonic  series  solution  involves  sums  and 
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differences  of  very  large  terms  since  the  constants  and  functions  are  pro¬ 
portional  to  powers  of  y.  ,  yo  ,  and  y.  This  situation,  of  course,  can  lead 
to  roundoff  errors  in  the  numerical  implementation;  the  appearance  of  these 
errors  depends  on  the  precision  of  the  machine  being  used.  Therefore,  the 
first  set  of  studies  which  was  carried  out  simply  determined  the  y-distance 
from  the  axis  of  symmetry  at  which  roundoff  errors  began  to  significantly 
affect  the  solution  when  evaluated  on  the  University  of  Illinois  CDC  Cyber 
175  digital  computer.  Both  single  precision  (14  digit  accuracy)  and  double 
precision  (29  digit  accuracy)  versions  of  the  subroutine  which  performs  all 
of  the  constant  and  function  evaluations  for  the  series  solution  were 
tested.  As  expected,  it  was  found  that  the  only  significant  parameter  in 
these  tests  was  y^  ,  the  dimensionless  distance  from  the  axis  of  symmetry  to 
the  inner  nozzle  wall.  This  quantity  may  be  used  to  characterize  the  axis- 
to-throat  distance  since  y  •-  y  <  y  and  y  -  y  =  1  in  the  region  of 
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interest.  Other  parameters  such  as  B,  the  angle  between  the  throat  axis 
and  the  symmetry  axis,  n,  the  parameter  in  the  expansion  variable  defini¬ 
tion,  R  and  R  ,  the  wall  radii  of  curvature,  etc.,  were  found  not  to 

c  i  co 

affect  tne  appearance  of  roundoff  errors. 

The  results  for  the  upper  limits  on  y.  for  avoiding  roundoff  errors  in 
the  determination  of  the  various  orders  of  perturbation  velocity  components, 
(ui  * vi  j  *  ^ u2  ’ v 2  j  ’  and  » v3 ] »  a|re  shown  in  part  (a)  of  Table  III. 2.  The 

single  precision  limits  were  found  by  noting  the  y  value  at  which  the 
single  precision  results  for  the  velocity  components  in  the  throat  plane 
deviate  significantly  [t>(0.1%)]  from  the  double  precision  values,  while  the 
double  precision  limits  were  found  by  observing  the  value  of  y.  at  which  the 
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aouble  precision  values  deviate  significantly  ~_(0.K)]  from  those  of  the 
planar  limit  which  is  asymptotically  approached  as  y  -  If  all  three 
orders  in  the  series  solution  are  employed,  the  limit  for  single  precision 
is  approximately  y.  =15;  this  includes  the  axi symmetric-no  centerbody  case 
and  most  annular  geometries  of  interest.  The  double  precision  limit  is 
y.  =  1000,  which  provides  a  very  good  approximation  to  the  plane  configura¬ 
tion.  If  only  two  terms  of  the  solution  are  used,  the  limits  are  y.  =  60 
and  y.  =  30,000,  respectively,  for  single  and  double  precision,  etc.  .  The 
reason  that  the  limits  are  relaxed  for  the  lower  orders  is  that  their  solu¬ 
tions  involve  lower  powers  of  y  than  does  the  third  order  solution. 

Since  the  derivation  of  the  discharge  coefficient,  C  ,  requires  the 
integration  of  the  density-u  velocity  component  product  across  the  throat 
plane,  the  constants  involved  in  its  evaluation  contain  higher  powers  of  y 
and  yo  than  do  the  velocity  components  themselves.  Thus,  the  upper  limits 
on  y.  for  avoiding  roundoff  errors  are  more  restrictive  for  CD  as  shown  in 
part  (b)  of  Table  III. 2.  The  technique  used  to  determine  these  limits  was 
similar  to  that  previously  discussed  for  the  perturbation  velocity  com¬ 
ponents  except  that  in  this  case  it  was  the  flow  coefficient  constants,  C  , 
CD2 ,  and  CD3  tnat  were  monitored.  For  solutions  involving  all  three  orders, 
the  limit  on  y.  is  approximately  10  for  single  precision  and  250  for  double 
precision.  However,  since  the  value  of  the  discharge  coefficient  is  so 
close  to  1.0  for  all  of  the  configurations  for  which  the  analysis  applies, 
i.e.,  for  Rc2l0(1):  .98  <_  CD  £  1.0,  the  discharge  coefficient,  with  its 

more  stringent  limits,  is  of  less  interest  than  the  velocity  components. 

This  is  particularly  true  since  one  of  the  primary  purposes  of  the  nozzle 
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throat  analysis  is  to  establish  an  accurate  initial  value  line  for  analyzing 
the  supersonic  portion  of  the  flowfield  using,  e.g.,  the  method  of 
character!'  sties . 

At  this  point  it  should  be  mentioned  that  the  numerical  implementation 
of  the  series  solution  is  extremely  fast  on  the  Cyber  175.  A  typical  prob¬ 
lem  requires  less  than  a  second  of  execution  time  when  the  single  precision 
version  uf  the  constant-evaluating  subroutine  is  used  and  less  than  two 
seconds  when  the  double  precision  version  is  employed.  For  this  reason  the 
double  precision  form  is  routinely  used,  and  as  can  be  deduced  from  the  pre¬ 
ceding  discussion,  it  provides  throat  flowfield  information  for  a  wide 
variety  of  configurations.  Conventional  axi symmetric,  annular,  and  planar 
geometries  can  all  be  analyzed  as  the  distance  from  the  axis  to  the  inner 
boundary  is  increased  from  zero  to  the  limiting  value  for  avoiding  roundoff 
errors. 

Once  these  limits  on  the  axis-to-throat  distance  were  established,  a 
second  series  of  studies  was  performed  to  observe  the  behavior  of  the 
transonic  series  solution  over  a  wide  range  of  parameters  for  nozzles  with 
circular  arc  contours.  The  parameters  involved  in  the  analysis  are:  y, 
the  specific  heat  ratio  of  the  flowing  gas;  n*  the  parameter  in  the  expansion 
variable  definition,  Eq.  (111-22);  the  geometrical  configuration  of  the 
nozzle  including  y.  ,  R  ,  Rc<},  and  8;  and  the  number  of  terms  of  the  expansion 
solution  to  be  included,  i.e.,  first,  second,  or  third  order  solutions. 

There  are  an  infinite  number  of  possible  geometries  and  values  of  y  and  n  as 
well  as  the  three  possible  solution  orders  which  can  be  investigated. 
Obviously,  then,  these  must  be  limited  in  a  significant  and  useful  way  in 
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order  to  effectively  study  the  solution  behavior.  First,  it  was  found  that 
the  solutions  are  not  strongly  dependent  on  y  so  that  only  the  usual 
diatomic  value  of  1.4  was  used.  Hence,  it  is  assumed  that  the  analysis 
correctly  predicts  the  weak  parametric  dependence  upon  this  physical 
property.  Five  values  of  n  were  employed,  n  =  0,  0.5,  1,  2,  and  4,  and 
each  of  the  three  solution  orders  was  considered.  In  addition,  twelve  total 
configurations  were  analyzed  with  three  values  of  the  radius  of  curvature, 

=  0.625,  2,  and  5,  utilized  in  each  of  four  groups  of  geometries.  These 
configurations  are  defined  and  shown  in  Fig.  III. 2:  (1)  conventional  axi- 
symmetric;  (2)  annular  axi symmetric  with  the  outer  boundary  straight  and 
parallel  to  the  axis  of  symmetry,  y.  =  2;  (3)  annular  with  the  inner  radius 
of  curvature  twice  that  of  the  outer  and  the  x-axis  inclined  at  roughly 
+  20°  to  the  axis  of  symmetry,  y.  =  5.3;  and  (4)  an  approximately  plane 
symmetric  configuration  with  y.  =  100.  The  value  R,  =  0.625  was  chosen 
since  a  dimensionless  wall  radius  of  curvature  less  than  unity  was  desired 
and  this  was  the  radius  of  curvature  value  for  the  axisyiimetric  nozzle 
tested  by  Cuffel ,  et  al.  [30].  As  mentioned  in  Chapter  II,  this  is  the  only 
experimental  investigation  to  date  in  which  flowfield  measurements  have  been 
obtained  for  an  axisymmetric  configuration. 

With  these  60  combinations  of  nozzle  geometry  and  n,  the  following 
studies  were  performed.  For  each  case,  the  M  =  0.8,  1.0,  and  1.2  constant 
Mach  number  contours  were  first  found  and  plotted  as  one,  two,  and  three 
terms  of  the  expansion  solution  were  included.  These  plots  provide  graphi¬ 
cal  information  concerning  the  convergence  properties  of  the  solution  over 
the  range  of  parameters  investigated.  The  convergence  was  further  studied 
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2  2  3  3 

by  printing  the  quantities  u, ,  v  -  ,  u2  -  <  v2  '  u3 '  *  and  v3  ■  at  a  number 
of  points  along  the  loci  of  the  third  order  iso-Mach  curves  and  noting 
whether  the  following  inequalities  were  satisfied:  i u, - '  u  , 
j  u3  e |  <  |u  !,  etc.  .  In  a  similar  manner,  the  discharge  coefficient  inequal¬ 
ities,  jC  t:  j  <  | C  j  and  |  C  c |  <  jC  j,  were  also  investigated  for  each 

L/  U I  Ls  a  LJ  £ 

configuration.  Satisfaction  of  these  inequalities  indicates  that  the  series 
solution  is  convergent  for  the  combination  of  parameters  under  consideration. 

The  final  set  of  tests  in  this  series  involved  evaluating  the  u  and  v 
velocity  components  from  the  first,  second,  and  third  order  solutions  at  a 
number  of  points  along  the  third  order  M  =  0.8,  1.0,  and  1.2  constant  Mach 
number  curves  and  substituting  these  values  into  the  exact  governing  equa¬ 
tions.  In  this  context  "exact"  is  used  to  denote  a  form  of  the  equations  in 
which  no  order  of  magnitude  simplifications  have  been  made  or  expansions 
assumed.  The  equations  which  were  utilized  in  this  study  are  irrotational ity 
condition  ( 1 1 1-8)  and  gas  dynamic  equation  ( 1 1 1- 9 ) ,  repeated  here  for 
convenience, 

u  -  v  =  0  ( II 1-8) 

y  x 


b2-&»2 


— 7T  UVU  + 
Y+l  y 


1-V 


Ill  u2 


y-1  2  Y-l  2 

1  -  u  -  77T  v 


Y+l 

(v  cosB  +  u  sing) 


=  0 


( 1 1 1 -9 ) 


y  cos  B  +  x  sing 

Whereas  the  identical  satisfaction  of  the  equations  and  boundary  conditions 
in  the  expansion  formulation,  Eqs.  (III-38)-(II  1-42 )  and  ( 1 1 1-48)- ( 1 1 1-53) , 
verifies  the  details  of  the  solution  development  as  discussed  in  the  preceding 
section,  the  size  of  the  residuals  in  the  exact  equations  above  provides  a 
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measure  of  the  validity  of  the  approximations  made  during  the  course  of  tne 
formulation.  Comparison  of  the  residuals  among  the  various  cases  studied 
also  allows  conclusions  to  be  drawn  concerning  the  effects  of  the  parameters 
on  the  various  assumptions  and  approximations  which  have  been  made.  As  in 
the  solution  verification  studies,  the  derivatives  in  the  exact  governing 
equations  have  been  approximated  by  second  order,  finite  differences  with  a 

.  -4 

grid  spacing  of  ix  =  Ay  =  10 

The  exact  boundary  conditions  that  the  velocity  vector  be  tangent  to 
the  bounding  walls  was  also  checked  along  each  contour  from  the  M  =  0.8  loca¬ 
tion  to  the  M  =  1.2  location.  This  interval  was  widened  in  order  to  include 
the  throat  station  for  those  cases  in  which  it  was  not  already  included. 

The  exact  boundary  conditions  were  investigated  in  the  form, 

=  0  ( 1 1 1-83) 

G(  R,  Z)  =  0 

=  0  ,  ( 1 1 1-84) 

H(  R,  Z)  =  0 

where  G(R,Z)  =  0  and  H(R,Z)  =  0  are  the  equations  of  the  inner  and  outer  con¬ 
tours  and  (0+6)  is  the  inclination  angle  of  the  velocity  vector  in  the 
cylindrical  R-Z  coordinate  system. 

The  basic  purpose  of  these  studies  was  to  investigate  in  an  organized 
and  rational  manner  the  convergence  properties  of  the  solution  and  to  deter¬ 
mine  if  an  optimum  value  of  the  parameter  n  could  be  found  that  provided 
improved  solution  convergence  and  simultaneously  minimized  the  residuals  in 
the  exact  governing  equations  and  boundary  conditions.  The  main  conclusions 
of  this  investigation  are  discussed  below. 


dR 

dZ 


tan (0+6) 


-  tan(0+6) 
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(1)  As  expected  from  the  problem  formulation,  all  orders  of  the  expan¬ 
sion  solution  identically  satisfy  both  the  exact  irrotational ity  condition 
and  the  boundary  conditions  for  straight  walls.  Therefore,  the  significant 
tests  are  the  residuals  for  the  exact  gas  dynamic  equation  and  for  the 
boundary  conditions  at  curved  contours,  as  well  as  the  convergence  prop¬ 
erties  of  the  series  solution. 

(2)  For  the  smallest  wall  radius  of  curvature  tested  in  each  of  the 
groups  of  geometries,  R^  =  0.625,  the  solutions  are  highly  divergent  for 

n  =  0,  Fig.  III. 3.  This  behavior  is  also  expected  since  for  n  =  0  the  solu¬ 
tions  correspond  to  those  of  Hall  [35]  who  used  the  expansion  parameter 
c  =  Rc  1  ,  and  for  Rc  <  1  this  expansion  parameter  exceeds  unity.  As  n  is 
increased  for  these  configurations,  the  convergence  of  the  solutions  is 
dramatically  improved,  Fig.  III. 4,  although  for  the  highest  value  tested 
here,  n  =  4,  the  solutions  tend  to  be  only  slowly  convergent,  especially  away 
from  the  sonic  line.  Interestingly,  as  n  is  increased  from  zero  for  these 
small  radius  of  curvature  geometries,  the  convergence  of  the  higher  Mach 
number  contours  is  improved  first. 

(3)  For  configurations  with  large  radius  of  curvature  boundaries,  the 
high  subsonic,  i.e.,  M  =  0.8,  constant  Mach  number  contours  are  less  con¬ 
vergent  than  the  low  supersonic,  i.e.,  M  =  1.2,  ones  for  n  =  0,  Fig.  III. 5. 

For  larger  values  of  n,  improved  convergence  of  the  high  subsonic  iso-Mach 
curves  is  obtained  somewhat  at  the  expense  of  the  convergence  of  the  low 
supersonic  ones.  Fig.  III. 6.  It  should  also  be  mentioned  that  for  nozzles 
with  a  large  throat  wall  radius  of  curvature,  there  is  little  difference 
between  the  third  order  constant  Mach  number  contours  for  various  values  of  n. 
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(4)  For  n  =  0,  the  expansion  solution  is  generally  an  alternating 
series,  such  that  the  third  order  iso-Mach  curve  for  a  particular  Mach  num¬ 
ber  lies  between  those  of  the  first  and  second  order.  Fig.  III. 7.  As  n  is 
increased,  the  nature  of  the  solution  gradually  changes  until  at  n  =  4  it  is 
a  slowly  convergent,  non-al ternating  type  series.  Fig.  III. 8. 

(5)  For  geometries  and  values  of  n  for  which  the  solutions  are  conver¬ 
gent,  the  residuals  in  the  gas  dynamic  equation  are  generally  reduced  as  the 
order  of  the  solution  is  increased  from  the  first  to  the  second  to  the  third 
order.  The  residuals  in  the  boundary  conditions  for  the  circular  arc 
boundaries  are  also  generally  reduced  as  the  higher  order  solutions  are 
utilized  with  this  trend  becoming  very  strong  as  n  is  increased  from  zero. 

(6)  It  appears  that  the  dependence  of  the  residuals  in  the  gas  dynamic 
equation  upon  n  is  directly  related  to  the  convergence  of  the  series  solu¬ 
tion.  Thus,  those  values  of  n  which  result  in  improved  convergence  prop¬ 
erties  for  the  expansion  solution  also  generally  result  in  reduced  residuals 
in  the  gas  dynamic  equation. 

(7)  The  residuals  in  the  exact  boundary  conditions  for  the  circular 
boundaries  are  increased  near  the  throat  location  as  n  is  increased  from 
zero,  particularly  for  the  larger  values  of  n.  This  behavior  is  consistent 
with  the  observations  from  Eqs.  (II 1-32 )  and  (III- 33 )  that  the  second 
derivatives  in  the  Maclaurin  series  expansions  for  the  wall  contours  are 
exactly  matched  only  for  n  =  0  and  that  the  approximations  deteriorate  pro¬ 
gressively  as  n  is  increased. 

(8)  As  found  by  previous  investigators,  the  geometrical  parameter 
which  is  by  far  the  most  influential  is  the  radius  of  curvature  of  the 
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bounding  walls.  For  nozzles  having  a  large  radius  of  curvature,  the  solution 
behavior  is  generally  improved  over  that  for  nozzles  with  sharp  throat 
curvature,  i.e.,  the  residuals  in  the  gas  dynamic  equation  and  boundary 
conditions  are  reduced  and  the  solutions  are  more  convergent.  This  behavior 
is,  of  course,  expected  since  the  expansion  parameter  has  been  defined  as 
c  =  |\+rij  *  •  The  parametric  effects  on  the  solution  behavior  of  the  distance 
of  the  nozzle  throat  from  the  axis  of  symmetry  and  the  inclination  of  the 
main  flow  x-direction  from  the  axis  are  generally  negligible  so  that  the 
dependence  of  the  solution  on  n,  the  solution  order,  and  the  wall  radius  of 
curvature  discussed  previously  are  repeated  in  each  of  the  four  groups  of 
geometries  sketched  in  Fig.  III. 2.  Two  exceptions  to  this  statement  are: 

(a)  For  cases  in  which  both  boundaries  are  curved  it 
appears  that  the  shift  in  solution  nature  from  an 
alternating  type  series  to  a  non-alternating  one 
occurs  at  lower  values  of  n  than  for  those  cases 
in  which  one  boundary  is  straight. 

(b)  For  the  smallest  radius  of  curvature  cases  tested  in 
each  of  the  four  geometrical  groups,  the  solution 
for  the  conventional  axi symmetric  configuration, 

y.  =0,  appears  to  be  the  best  behaved.  For  the 
other  three  sharp  curvature  cases,  the  M  =  0.8,  and 
in  two  cases,  the  M  =  1.0  contours  were  found  not  to 
intersect  the  highly  curved  wall,  Fig.  III. 9. 

Based  on  the  preceding  observations,  it  is  concluded  that  third  order 
solutions  should  generally  be  used.  Also,  for  configurations  with  a  large 
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wall  radius  of  curvature,  the  expansion  solution  is  essentially  independent 
of  n.  However,  for  geometries  with  small  radius  of  curvature  walls,  two 
competing  effects  appear:  as  n  is  increased  from  zero  the  convergence  of 
the  series  solution  is  improved  significantly,  but  the  satisfaction  of  the 
exact  boundary  conditions  in  the  throat  region  is  compromised,  especially 
for  the  largest  value  tested,  n  =  4.  Ideally,  then,  the  smallest  value  of  n 
for  which  the  convergence  properties  in  the  region  of  interest  are  deemed 
satisfactory  should  be  employed.  For  general  use  n  =  2  is  recommended/ 

Third  order,  n  =  2  solutions  are  compared  to  two  existing  data  sets  in 
Figs.  III. 10  and  III. 11.  In  the  first  figure  comparison  is  made  to  the  data 
of  Cuffel ,  et  al.  [30]  for  an  axisymmetric  nozzle  with  a  dimensionless  wall 
radius  of  curvature  of  Rc  =  0.625.  As  mentioned  in  Chapter  II,  this  data 
was  obtained  by  means  of  a  small  diameter,  hollow  tube  which  was  stretched 
between  an  upstream  stagnation  chamber  and  a  downstream  vacuum  chamber. 

This  small  value  or  Rc  provides  a  severe  test  for  the  expansion  solution,  and 
as  can  be  seen  in  the  figure,  the  results  are  in  reasonably  good  agreement 
with  the  data  through  the  throat  region  except  for  the  high  subsonic  Mach 
numbers,  M  =  0.6  and  M  =  0.8.  In  particular,  the  M  =  1.2  contour  appears  to 
provide  an  accurate  starting  line  for  initiating  the  hyperbolic  computations 
for  the  supersonic  region.  In  Fig.  III. 11  the  results  of  the  series  solution 
are  compared  to  the  sonic  line  data  of  Flack  and  Thompson  [14]  for  a  planar, 
symmetric  nozzle  with  a  wall  curvature  of  Rc  =  2.  The  experimental  informa¬ 
tion  was  determined  by  using  static  pressure  measurements  obtained  from 

tThe  equivalent  recommended  value  for  the  simplified,  conventional  axisym¬ 
metric  solution,  Appendix  C,  is  n  =  1. 
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pressure  taps  located  in  the  sidewalls  of  the  nozzle.  For  the  theoretical 
results,  the  plane,  two-dimensional  geometry  has  been  approximated  by  an 
annular  configuration  located  a  very  large  dimensionless  distance  from  the 
axis  of  symmetry,  y.  =  1000.  As  shown  in  the  figure,  the  analytical  and 
experimental  sonic  lines  are  of  the  same  shape  with  the  measurements  lying 
somewhat  downstream  of  the  series  results.  The  agreement  is  within  approxi¬ 
mately  3  percent,  however,  since  calculation  of  the  Mach  numbers  along  the 
experimental  sonic  line  using  the  expansion  solution  results  in  values  of 
about  1.03. 

As  a  result  of  the  parametric  studies  and  comparisons  presented  in  this 
section,  it  is  felt  that  the  approximate  solution  which  has  been  developed 
is  applicable  to  a  wide  variety  of  nozzle  configurations  including  axisym- 
metric,  annular,  and  planar  ones.  As  long  as  the  radius  of  curvature  of  the 
bounding  walls  is  of  the  order  of  the  throat  half-height  or  larger,  and  as 
long  as  attention  is  restricted  to  the  transonic  throat  region,  x  =  o(e1//2), 
with  the  inclination  angle  satisfying  the  restriction,  B  /y  =  (7(1 ),  satis¬ 
factory  results  are  expected.  Further  comparisons  between  the  theoretical 
solution  and  experimental  measurements  obtained  during  the  course  of  this 
investigation  are  presented  in  the  next  chapter. 
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Figure  III.l  Configuration  for  throat  flowfield  analysis  of 
annular  supersonic  nozzles 
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Table  III.l  Solution  Tests 


I.  Reduction  to  Previous  Solutions  for  Simpler  Geometries 


A.  Reduction  to  Conventional  Axisynmetri c  Configuration  as  y.  -  0 

1.  n  =  0  --  Hall  [35],  with  corrections  by  Kliegel  and  Levine  [46]. 

2.  n  =  2  —  Dutton  [87] . 

B.  Reduction  to  Planar  Symmetric  Configuration  as  y,  -*■  ^ 

1.  n  =  0  --  Hall  [35] . 

2.  n  =  2  --  Dutton  [87]  and  Thompson  and  Flack  [37], 

C.  Reduction  to  Annular  Ax i symmetric  Configuration  with  One  Boundary 
Straight  and  Parallel  to  the  Axis  of  Symmetry. 

1.  n  =  2,  g"(0)  =  0  --  Dutton  [87]. 

2.  n  =  2,  h" (0)  =  0  --  Dutton  [87]. 


II.  Independent  Numerical  Back  Substitution  into  the  Governing  Equations  and 
Boundary  Conditions  (Derivatives  approximated  with  finite  differences) 
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Table  III. 2  Roundoff  Error  Investigation 


(a)  Approximate  values  of  y  at  which  roundoff  error*  affects  solutions 
for  perturbation  velocity  components: 


Single  Precision 

Double  Precision 

First  Oraer 

h-.i 

3000 

io7 

Second  Oraer 

60 

30,000 

Third  Order 

h'vd 

15 

1000 

(b)  Approximate  values  of  y  at  which  roundoff  error*  affects  solutions 
for  discharge  coefficient  constants: 


Single  Precision 

Double  Precision 

First  Order 

M 

120 

50,000 

Second  Order 

M 

20 

1500 

Third  Order 

(CD3  ) 

10 

250 

*0n  the  University  of  Illinois  CDC  Cyber  175  digital  computer 


(a)  Conventional  axisynmetric 


Figure  III. 2  Configurations  for  parametric  study  of  solution  behavi 


(d)  Approximately  plane  symmetric,  y.  =  100 

Figure  III. 2  (cont.)  Configurations  for  parametric  study  of 

solution  behavior 


Figure  111.3  Sonic  lines  for  axisymmetric  configuration 
R  =  0.625,  n  =  0 


I 

I 
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utr s  f nr  annular  configuration. 
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Figure  III. 6  M  =  0.8,  1.0,  and  1.2  contours  for  annular  configuration, 
Rri  =  5’  =  2,  n  =  2 
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Figure  III. 7  M  =  0.8,  1.0,  and  1.2  contours  for  approximately  planar 
configuration,  R  =  2,  y.  =  100,  n  =  0 

CO  1 


MACH 

NUMBER  CONTOURS 

7  =1.4 

Y,=0.0 

71  =2.0 

Rcr  00 

0  =0.0 

Roq=.625 

SYMBOLS  =  DATA  OF  CUFFEL,  ET.  AL.C30] 

-- 

—  SERIES  SOLUTION 

Figure  III. 10  Comparison  of  third  order,  n  =  2  solutions  to  data 
of  Cuffel,  et  al.  [30]  for  an  axi|ymmetric  nozzle 
with  R  =  0.625;  ReoJ  =  2.8  x  10  for  experiments 
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Figure  III. 11  Comparison  of  third  order,  n  =  2  solutions  to  data 
of  Flack  and  Thompson  [14]  for  a  planar  nozzle  with 
R  =  2.0;  average  Re2d  =  8.6  x  10  for  experiments 
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IV.  EXPERIMENTAL  INVESTIGATION 

As  discussed  in  the  literature  survey,  there  have  been  few  experimental 
studies  of  supersonic  nozzle  throat  flows,  particularly  in  terms  of  flow- 
field  measurements  for  axisymmetric  configurations.  For  this  reason,  an 
experimental  investigation  was  performed  to  obtain  flowfield  pressure  dis¬ 
tributions  for  a  number  of  annular  nozzle  geometries  including  a  conventional 
axisymmetric  nozzle.  For  the  annular  configurations,  these  measurements  are 
the  first  ones  known  to  consider  more  than  just  the  wall  pressure  distribu¬ 
tion  along  one  of  the  bounding  walls,  while  the  data  of  Cuffel ,  et  al .  [30] 
is  the  only  other  set  of  flowfield  data  for  the  axisymmetric,  no-centerbody 
case.  The  resulting  measurements  provide  a  basis  of  comparison  in  order  both 
to  verify  the  theoretical  solution  developed  in  the  preceding  chapter  and  to 
draw  conclusions  regarding  its  range  of  applicability. 

A.  EQUIPMENT  AND  APPARATUS 

The  main  objections  to  the  stretched  tube  technique  employed  in  [30]  to 
obtain  the  nozzle  throat  flowfield  pressure  distributions  are  that  the  tube 
introduces  disturbances  of  unknown  magnitude  into  the  flowfield  and  for  off- 
axis  positions  the  velocity  vector  is  not  necessarily  parallel  to  the  tube. 
Rather  than  using  this  method,  the  one  which  has  been  utilized  in  this  study 
is  the  splitter-plate  technique  whereby  half-section  cylindrical  models  were 
constructed  and  mounted  on  a  plate  whose  surface  corresponds  to  a  symmetry 
plane  for  the  axisymmetric  geometry.  Pressure  taps  were  placed  on  the 
splitter  plate  in  the  transonic  region  of  interest  from  which  the  desired 
transverse  and  streamwise  pressure  distributions  were  obtained. 
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The  splitter  plate  also  introduces  disturbances  into  the  flowfield  in 
the  form  of  the  boundary  layers  which  grow  both  on  the  plate  and  in  the 
corners  where  the  half-models  meet  the  plate.  These  boundary  layers,  of 
course,  are  not  present  for  the  full  axisyrmetri c  configuration.  However, 
for  the  high  Reynolds  number,  favorable  pressure  gradient,  nozzle  flows  under 
consideration,  the  boundary  layers  in  the  throat  region  are  extremely  thin. 

In  fact,  sample  calculations  using  the  code  developed  in  [88]  show  that  for 
the  models  and  typical  operating  conditions  used  here,  the  boundary  layer 
displacement  thickness  in  the  throat  region  is  less  than  0.05  mm  which  is 
less  than  0.2%  of  the  throat  diameter  in  the  worst  case.  It  is  therefore 
felt  that  the  measurements  obtained  with  the  half-section  models  correspond 
closely  to  the  "true"  measurements  for  the  full  axisymmetric  geometries.  To 
further  test  this  hypothesis  both  a  full-section  and  the  corresponding  half¬ 
section  models  were  constructed  for  a  conventional  axisymmetric  nozzle  con¬ 
figuration.  Comparison  of  the  wall  pressure  measurements  from  the  two  models 
then  allows  conclusions  to  be  drawn  regarding  the  effects  of  the  splitter 
boundary  layers.  The  results  of  these  experiments  will  be  reported  in  a 
later  section. 

Figure  IV. 1  is  a  photograph  of  the  full  section  axisymmetric  nozzle 
which  was  designed  and  constructed  showing  the  wall  pressure  taps  which  were 
spiraled  around  the  circumference  of  the  nozzle  through  the  throat  region. 
Also  visible  in  this  photograph  are  the  last  few  static  holes  which  were 
carefully  drilled  normal  to  the  wall  contour;  all  burrs  and  irregularities 
were  removed  at  the  wall  in  the  vicinity  of  the  pressure  taps.  Unless 
specifically  stated  otherwise,  the  static  holes  all  have  the  same  diameter 
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of  0.51  mm  and  the  model  material  is  aluminum.  Figure  IV. 2  shows  the  corre¬ 
sponding  half-section  axisymmetric  model.  The  half-nozzle  in  the  right 
background  is  mounted  on  the  splitter  plate  in  the  foreground  where  the  flow 
direction  is  from  left  to  right.  The  entrance  plate,  which  forms  the  left 
half  of  the  splitter  plate,  fits  inside  an  existing  axisymmetric  converging 
nozzle  whose  elliptic  entrance  section  matches  the  curved  boundaries  of  the 
entrance  plate.  The  blocking  flange  in  the  left  background  is  attached  and 
sealed  to  the  bottom  of  the  splitter  plate  so  that  no  flow  is  allowed  to 
leak  through  the  test  section  under  the  plate.  In  addition  to  the  8  wall 
taps  on  the  half-nozzle  which  were  used  to  obtain  the  wall  pressure  distribu¬ 
tion  for  comparison  with  that  of  the  full  nozzle,  71  static  pressure  taps 
were  located  on  the  splitter  plate  so  that  the  flowfield  pressure  distribu¬ 
tions  could  be  measured  for  this  axisymmetric  nozzle.  These  taps  have  been 
arranged  in  12  staggered  rows  thereby  providing  throat  region  pressure 
measurements  in  a  square  grid,  3.81  mm  on  a  side.  It  should  also  be  mentioned 
that  when  the  experiments  were  performed,  the  screw  heads  seen  on  the  surface 
of  the  entrance  plate  were  filled  in  with  epoxy  and  sanded  so  that  the  entire 
surface  of  the  splitter  plate  was  extremely  smooth  and  highly  polished. 

In  order  to  obtain  data  for  annular  configurations,  the  centerbody  shown 
in  Fig.  IV. 3  is  added.  It  should  be  noted  that  in  terms  of  the  theory  pre¬ 
sented  in  the  last  chapter,  only  the  central  portion  of  the  centerbody,  whose 
radius  of  curvature  in  the  meridional  plane  is  50.8  mm,  is  of  importance. 

The  cylindrical  sections  upstream  and  downstream  of  this  region  merely  pro¬ 
vide  a  means  for  attaching  the  centerbody  to  the  splitter  plate.  By  displac¬ 
ing  the  centerbody  relative  to  the  outer  nozzle,  measurements  for  various 


84 


annular  configurations  could  be  made.  As  will  be  discussed  in  more  detail  in 
a  later  section,  three  such  configurations  were  investigated.  The  splitter 
plate  in  this  case  is  not  the  same  one  used  for  the  conventional  axisymmetric 

geometry  of  Fig.  IV. 2  since  a  slot  is  required  for  the  centerbody  wall  taps 

and  dowel  pin  and  bolt  holes  are  required  to  locate  and  hold  the  centerbody 
on  the  plate.  The  pressure  taps  on  the  splitter  plate  in  this  case  have 
been  arranged  in  10  alternating  rows  of  5  taps  each  so  that  pressure  measure¬ 
ments  in  a  square,  3.18  mm  grid  are  obtained.  The  measurements  from  these 
50  splitter  plate  taps  together  with  the  9  centerbody  wall  taps  and  the  8 
nozzle  wall  taps  allow  the  pressure  field  in  the  annular  nozzle  throat  to  be 
completely  mapped. 

Two  views  of  the  assembled  apparatus  of  Fig.  IV. 3  are  shown  in  Figs.  IV. 4 
and  IV. 5.  In  Fig.  IV. 4,  the  view  is  in  the  downstream  flow  direction  through 
the  test  section.  In  viewing  this  photograph,  it  should  be  remembered  that 

the  entrance  plate  is  highly  polished  aluminum  so  that  everything  seen  on  its 

surface  is  a  reflection  from  the  half-nozzle  mounted  on  the  splitter  plate. 

The  pressure  tap  grid  is  detailed  in  Fig.  IV. 5,  which  looks  upstream  through 
the  measurement  region.  The  shorter  row  of  taps  seen  below  the  plate  is  the 
set  which  passes  through  the  slot  in  the  splitter  plate  and  is  attached  to 
the  centerbody  for  measurement  of  its  wall  pressure  distribution.  In  assembl¬ 
ing  the  half-section  models,  Loctite  anaerobic  adhesive  was  used  on  all  of  the 
flat  mating  surfaces  and  silicone  rubber  sealant  was  also  employed  along  all 
seams  in  noncritical  regions.  Together  with  the  0-rings  which  were  utilized 
between  mating  flanges,  this  sealing  method  prevented  the  occurrence  of  leaks 
in  the  half-model  test  sections. 
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A  photograph  of  the  partially  assembled,  annular,  supersonic  nozzle 
test  section,  as  it  is  mounted  on  the  brass  converging  inlet  nozzle,  is  given 
in  Fig.  IV. 6.  The  converging  nozzle,  in  turn,  is  bolted  to  the  faceplate  of 
one  of  the  test  stands  in  the  supersonic  wind  tunnel  area  of  the  Mechanical 
Engineering  Laboratory.  However,  t he  supersonic  jet  cannot  simply  be 
exhausted  to  ambient  conditions  because  of  the  very  high  noise  levels 
associated  with  such  a  jet.  Instead,  the  test  section  is  connected  to  the 
facility  silencing  system.  In  this  arrangement  the  supersonic  jet  exhausts 
into  a  viewing  chamber  which  also  contains  a  conically  convergent  jet- 
catcher/diffuser.  A  view  of  the  annular  nozzle  test  apparatus  as  seen 
through  one  of  the  side  windows  of  the  exit  chamber  is  shown  in  Fig.  IV. 7. 
Downstream  of  this  chamber  the  jet  flow  enters  a  double-walled  silencing 
duct  before  passing  through  a  final  silencer  and  exiting  to  atmosphere  out¬ 
side  of  the  laboratory.  An  overall  view  of  the  assembled  test  apparatus  is 
given  in  Fig.  IV. 8  showing  the  first  few  feet  of  the  double-walled  silencing 
duct  at  the  right.  Dry,  filtered,  compressed  air  is  supplied  to  the  stagna- 

3 

tion  test  chamber  by  a  Gardner-Denver  screw  compressor  through  a  140  m  tank 
farm  and  a  152.4  mm  supply  line.  The  pressure  in  the  stagnation  chamber  is 
regulated  by  means  of  a  Fisher-Governor  automatic  controller  and  control 
valve.  A  schematic  of  the  air  flow  circuit  is  presented  in  Fig.  IV. 9. 

Pressures  were  measured  with  two  CEC  type  4-312  unbonded  strain  gage 
pressure  transducers.  A  0-700  kPa  gage  transducer  was  utilized  in  conjunc¬ 
tion  with  a  48-port,  motor-driven  model  J  Scanivalve  to  record  the  test  sec¬ 
tion  static  pressures.  The  stagnation  chamber  pressure  was  also  continuously 
monitored  by  means  of  a  0-1000  kPa  gage  transducer.  These  pressure  signals 
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were  conditioned  with  two  Daytronic  model  870  conditioner-amplifier  modules. 
The  resulting  output  signals  were  recorded  from  a  Oaytronic  model  890  digi¬ 
tal  meter  and  a  Keithley  model  179  digital  voltmeter  for  the  stagnation  and 
static  pressures,  respectively.  Additional  equipment  used  in  the  course  of 
the  experiments  is  also  described  in  the  next  section. 

B.  PROCEDURE 

The  first  step  in  the  experimental  procedure  involved  calibration  of  the 
pressure  transducers.  This  was  accomplished  by  using  vacuum  and  pressure 
sources  in  conjunction  with  two  pressure  measurement  standards:  a  Wallace 
and  Tiernan  precision  mercury  manometer  for  the  0-200  kPa  absolute  pressure 
range  and  a  Wallace  and  Tiernan  bourdon  tube  absolute  pressure  gage  for  the 
200-500  kPa  range.  Of  the  two  devices  the  precision  manometer  is  the  more 
accurate,  but  it  is  l.mited  to  the  lower  pressure  range  by  the  allowable 
height  of  the  mercury  column.  The  calibrations  were  completed  by  obtaining 
linear,  least-squares  regressions  between  the  digital  readings  of  the  output 
devices  and  the  readings  from  the  precision  manometer  and  the  absolute  gage. 
For  both  transducers  the  calibration  data  was  found  to  be  extremely  linear. 

After  the  hardware  for  a  given  model  configuration  had  been  assembled, 
the  experimental  procedure  for  each  test  consisted  first  of  recording  the 
atmospheric  pressure  from  a  Setra  Systems  electronic  barometer.  The  stagna¬ 
tion  supply  pressure  was  then  set  with  the  Fisher-Governor  automatic  con¬ 
troller  and  the  static  and  stagnation  chamber  pressures  read  from  the  digital 
output  devices  as  the  Scanivalve  was  stepped  from  port  to  port  with  its 
solenoid  controller.  Since  the  half-section  models  all  had  more  than  48 
static  taps,  which  is  the  number  of  available  Scanivalve  ports,  approximately 


half  of  the  taps  were  connected  to  the  Scanivalve  tor  a  given  run  and  the 
rest  were  blocked  off.  In  later  tests  these  connections  were  reversed.  The 
supply  pressure  in  all  of  the  experiments  was  approximately  400  kPa.  It  was 
not  possible  to  vary  this  pressure  level  over  a  wide  range  since  it  had  to 
be  high  enough  to  ensure  that  the  measurement  region  was  free  from  the  shock 
waves  which  occur  due  to  the  exit  boundary  conditions  but  low  enough  so  that 
sufficient  time  i<is  available  to  record  the  data  in  these  blowdown  tests. 

For  these  reasons  the  Reynolds  number  based  on  sonic  conditions  and  the 
throat  height,  Re2d  ,  also  was  not  varied  over  a  wide  range. 

To  check  tie  repeatability  of  the  experiments,  five  runs  were  made  for 
each  geometrical  arrangement.  For  each  tap  the  data  was  reduced  by  calculat¬ 
ing  the  mean  and  standard  deviation  of  the  stati c-to-stagnation  pressure 
ratio  readings,  p/p0 ,  for  the  five  tests.  The  resulting  mean  values  were 
used  in  all  of  the  data  presented  in  the  next  section.  In  all  cases  the 
experiments  were  found  to  be  highly  repeatable.  The  standard  deviation  for 

_  3 

a  given  measurenent  was  generally  less  than  10  in  units  of  p/p0  where  this 
variable  spans  almost  the  entire  range  from  unity  to  zero  through  the 
transonic  throat  region.  This  corresponds  to  a  Mach  number  standard  devia- 

_  3 

tion  of  1.6  x  10  at  sonic  conditions. 

Before  presenting  the  experimental  results  it  is  appropriate  to  discuss 
the  accuracy  of  the  equipment  used  to  obtain  the  pressure  measurements.  By 
utilizing  micrometers  and  other  precision  measurement  gages,  the  critical 
dimensions  of  the  air  flow  models,  including  the  static  hole  locations,  were 
checked  and  found  to  be  within  +  0.05  mm  of  the  design  specifications.  Since 
the  throat  diameters  of  the  models  range  from  approximately  32  to  50  mm  and 
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since  the  static  holes  are  0.51  mm  in  diameter,  the  inaccuracies  incurred  in 
the  precision  machinework  may  be  considered  negligible.  For  the  pressure 
measurements,  the  least  accurate  device  used  in  the  calibration  procedure 
and  the  determination  of  atmospheric  pressure  was  the  0-1000  kPa  bourdon 
tube  absolute  pressure  gage.  The  accuracy  of  this  device  is  listed  as 
+_  0.7  kPa.  Considering  the  additional  effects  of  calibration  nonlinearity, 
amplifier  drift,  measurement  fluctuations,  etc.,  it  is  felt  that  the  com¬ 
bined  accuracy  of  the  pressure  measurement  equipment  is  within  +  1.5  kPa. 

However,  since  the  data  is  generally  presented  in  terms  of  the  Mach 
number  rather  than  the  pressure  ratio,  it  is  natural  to  inquire  how  the 
inaccuracies  in  the  individual  measurements  of  the  static  and  stagnation 
pressures  affect  the  determination  of  the  Mach  number.  This  question  can  be 
answered  by  utilizing  functional  uncertainty  analysis  which  is  a  technique 
whereby  the  fractional  uncertainty  of  a  functional  relationship  is  calculated 
from  the  uncertainties  of  the  individual  measurements  involved  in  it.  As 
such  it  excludes  the  random  errors  of  particular  measurements  and  is  concerned 
only  with  the  accuracy  of  the  measuring  devices,  or  systematic  error.  The 
fractional  uncertainty  of  a  function  F  which  depends  on  i  variables,  F  =  F^x.  j, 
is  given  by  the  expression  [89], 

U(F)  _  _  fjZnL  UfXi  1 
F  i  Unxi  x. 

4 

where  U  is  used  to  denote  the  uncertainty.  For  the  case  at  hand,  the  isen- 
tropic  relation. 
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is  employed  to  convert  the  pressure  ratio  to  the  corresponding  Mach  number. 
Using  logarithmic  differentiation,  it  is  easily  shown  that, 


=  _  3 1 >:M  =  1  *  ^ 

a  Ik  p„  3’>:p  ,,2 

po  K  -vM 


( I V- 3 ) 


Therefore,  the  fractional  uncertainty  for  the  Mach  number  may  be  determined 
from. 


Using  1.5  kPa  as  the  uncertainty  of  both  the  static  and  stagnation  pressure 
measurements  and  400  kPa  as  the  value  of  the  stagnation  pressure  for  these 
experiments,  the  corresponding  fractional  uncertainty  of  the  Mach  number  can 
be  calculated  through  the  transonic  region.  For  M  =  0.6  a  value  of  1.3%  is 
obtained,  while  for  M  =  1.4  the  fractional  uncertainty  in  the  Mach  number  is 
0.63%.  Similarly,  the  fractional  uncertainty  for  the  static- to-stagnation 
pressure  ratio,  p/pQ ,  is  found  to  be  0.61%  at  M  =  0.6  and  1.25%  at  M  =  1.4. 
Thus  the  inaccuracies  in  the  determination  of  both  the  Mach  number  and  pres¬ 
sure  ratio  attributable  to  the  measuring  devices  are  expected  to  be  on  the 
order  of  one  percent  for  the  flowfields  investigated  here. 


C.  COMPARISON  OF  EXPERIMENTAL  AND  THEORETICAL  RESULTS 

Experimental  measurements  for  a  conventional  axisymmetric  configuration 
and  three  annular  ones  were  obtained.  A  half-section  drawing  of  the  particu¬ 
lar  annular  arrangement  for  which  the  centers  of  curvature  of  both  the  inner 
and  outer  wall  contours  lie  along  the  Z  =  0  plane  is  presented  in  Fig.  IV. 10. 
All  dimensions  are  in  millimeters.  The  air  flow  enters  the  test  section  from 
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the  stagnation  chamber  at  the  left  where  it  is  first  accelerated  through  the 
converging  nozzle  with  its  elliptical  entrance  section.  At  the  exit  of  the 
converging  nozzle  the  flow  cross-sectional  area  has  been  reduced  to  a 
50.8  mm  radius.  Further  acceleration  then  occurs  in  the  annular  supersonic 
nozzle.  Along  the  outer  nozzle  wall,  a  constant  area  cylindrical  region  is 
followed  by  37.5°  conically  convergent  section  which  is  tangent  to  a  circular 
arc  contour.  In  the  meridional  plane,  the  radius  of  curvature  of  this  arc 
is  25.4  mm,  and  the  radial  distance  to  the  minimum  point  on  the  outer  wall 
contour  is  also  25.4  mm.  Along  the  centerbody,  the  flow  is  first  accelerated 
from  the  stagnation  point  through  a  25.4  mm  radius  arc  to  a  cylindrical 
region  with  a  6.4  mm  radius.  This  section  is  followed  by  a  circular  arc  con¬ 
tour  whose  radius  of  curvature  is  50.8  mm  and  a  second,  downstream  cylindri¬ 
cal  region  with  a  6.4  mm  radius.  The  radial  distance  to  the  maximum  point  on 
the  centerbody  contour  is  9.5  mm.  For  the  axi symmetric,  Laval  nozzle  case, 
the  dimensions  for  the  outer  nozzle  are  identical,  and  the  centerbody  is 
simply  removed.  The  other  two  annular  configurations  are  obtained  by  dis¬ 
placing  the  centerbody  relative  to  the  outer  nozzle  both  upstream  and  down¬ 
stream  by  12.7  mm.  For  these  two  configurations  the  main  flow  x-direction  is 
inclined  with  respect  to  the  Z-axis  of  symmetry. 

To  facilitate  comparisons  between  the  various  configurations,  the  coordi¬ 
nate  system  used  in  the  presentation  of  the  data  is  the  R-Z  cylindrical  system 
shown  in  Fig.  IV. 10.  In  this  case  the  radial  R-axis  passes  through  the  mini¬ 
mum  point  on  the  outer  wall  contour  and  the  non-dimensionalizing  length  is 
chosen  as  the  radial  distance  to  this  minimum  point.  This  distance  is  also 
the  same  as  the  radius  of  curvature  of  the  outer  wall,  25.4  mm  (or  1  inch). 
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In  Fig.  IV. 11  the  wall  measurements  from  the  full-  and  half-section 
models  for  the  conventional  axisymmetric  configuration  are  compared.  The 
pressure  measurements  were  converted  to  Mach  numbers  using  the  isentropic 
relation,  Eq.  (IV-2),  and  for  this  case,  the  wall  radius  of  curvature  is 
equal  to  the  throat  half-height,  R  =1.0.  As  can  be  seen  in  the  figure, 

C 

the  wall  Mach  number  measurements  for  the  two  models  are  essentially  identi¬ 
cal,  with  the  root-mean-square  difference  between  the  data  sets  calculated  to 

—  2 

be  (AM)  =  0.89  x  10  .  For  the  first  two  and  last  four  pressure  taps  in 

r  ms 

the  streamwise  direction,  the  Mach  numbers  measured  with  the  half-section 
model  are  slightly  higher,  while  at  the  remaining  two  axial  stations, 

Z  =  -  0.3  and  Z  =  -  0.1b,  the  Mach  number  data  for  the  full  axisymnetric 
model  is  higher.  These  small,  somewhat  random  differences  between  the  mea¬ 
surements  could,  in  part,  be  due  to  slight  differences  in  the  pressure  tap 
locations  for  the  two  models  since  precise  location  of  the  taps  on  the  curved 
walls  was  difficult  both  to  achieve  and  to  check.  In  any  event,  the  very 
near  coincidence  of  the  two  data  sets  lends  credence  to  the  argument  that  the 
measurements  from  the  half-models  should  closely  approximate  those  for  the 
corresponding  full  axisymmetric  geometries.  Due  to  its  thinness,  the  bound¬ 
ary  layer  on  the  splitter  plate  apparently  has  little  effect  on  these 
measurements. 

Also  shown  in  Fig.  IV. 11  are  the  inviscid,  wall  Mach  number  results 
obtained  from  one-dimensional,  isentropic,  flow  with  area  change  analysis 
and  also  from  the  series  solution  developed  in  Chapter  III.  As  in  all  of  the 
comparisons  to  be  made  in  this  section,  the  third  order,  ri  =  2  expansion 
results  are  utilized  as  recommended  in  section  III.D.  From  somewhat  upstream 
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of  the  throat  location,  Z  =  0,  on  downstream,  it  is  seen  that  the  one¬ 
dimensional  analysis  is  seriously  in  error  while  the  series  solution  results 
agree  fairly  well  with  the  measurements.  However,  at  the  two  most  upstream 
wall  tap  locations,  Z  =  -  0.6  and  Z  =  -  0.45,  the  one-dimensional  results 
agree  more  closely  with  the  data  th  n  do  the  results  of  the  expansion  solu¬ 
tion.  This  lack  of  agreement  between  the  series  solution  and  the  wall  data 
at  these  upstream  stations  can  be  attributed  to  two  approximations  made  in 
the  theoretical  analysis.  First,  since  the  transonic  approximations  embodied 
in  Eqs.  (III-10)  and  (III-ll)  have  been  utilized,  it  is  expected  that  the 
solution  will  be  valid  only  for  nearly  sonic  conditions,  i.e.,  M  =  1. 
Therefore,  for  the  relatively  low  Mach  numbers  measured  at  the  two  upstream 
taps,  M  =  0.5,  it  is  not  surprising  that  the  series  solution  results  do  not 
agree  well  with  the  data.  On  the  supersonic  side,  however,  it  is  seen  that 
the  expansion  solution  agrees  with  the  wall  measurements  for  Mach  numbers 
well  above  1.5.  Probably  more  important  is  the  assumption  that  the  coeffi¬ 
cients  in  the  Maclaurin  series  expansion  for  the  wall  contours  have  the  same 
order  of  magnitude  estimates  as  those  for  circular  arcs.  From  the  axial 
location  Z  =  -  0.6  downstream  to  the  exit,  the  outer  wall  contour  of  the 
experimental  model  is,  in  fact,  a  circular  arc  with  R  =1.0.  However, 

C 

upstream  from  this  station,  the  wall  contour  of  the  model  consists  of  a 
37.5°  conically  convergent  section.  Hence,  it  would  not  be  expected  that 
wall  measurements  taken  near  the  tangency  point  between  the  conical  inlet 
and  the  circular  arc  contour  would  agree  well  with  the  theory.  The  transi¬ 
tion  in  the  contour  from  a  conical  to  a  circular  section  is  not  allowed  for 
in  the  analysis  since  no  provision  is  made  for  the  discontinuity  in  the 
second  derivative  which  occurs  at  the  tangency  point. 
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Figure  IV. 12  presents  a  comparison  between  the  contours  of  constant 
Mach  number  obtained  from  the  expansion  solution  and  those  determined  from 
the  splitter  plate  measurements  of  the  half-section  model  for  the  axisym- 
metric,  Laval  nozzle  under  consideration,  R  =1.0.  The  data  points  were 

C 

found  by  first  converting  the  stagnation  supply  pressure  measurement  and  the 
individual  static  tap  pressure  readings  to  the  corresponding  local  Mach  num¬ 
bers  using  Eq.  (IV-2).  Since  the  constant  Mach  number  contours  for  the 
configurations  studied  here  are  primarily  radial  in  nature,  the  data  points 
were  determined  by  finding  the  axial  Z-location  at  which  tht  given  Mach  num¬ 
ber  was  attained  along  each  of  the  pressure  tap  rows  that  is  a  constant 
radial  distance  from  the  axis  of  symmetry.  In  addition  to  the  outer  wall, 
there  are  seven  such  rows:  R  =  0,  0.15,  0.3,  0.45,  0.6,  0.75,  and  0.9. 
Hermite  spline  interpolation  was  utilized  in  obtaining  the  loci  of  the  points 
at  which  the  Mach  number  values  of  interest  occur  from  the  measurements  of 
the  individual  taps  along  each  row. 

As  can  be  seen  in  the  figure,  the  correspondence  between  the  series 
solution  results  and  the  experimental  measurements  is  generally  quite  good 
except  for  the  M  =  0.6  contour  and  near  the  wall  for  the  M  =  0.8  iso-Mach 
curve.  The  disagreement  in  these  regions  can  again  be  attributed  to  the 
fact  that  upstream  of  the  Z  =  -  0.6  location  the  wall  contour  of  the  experi¬ 
mental  model  is  conically  convergent  rather  than  being  the  circular  arc 
contour  of  the  downstream  region  as  used  in  the  implementation  of  the  expan¬ 
sion  solution.  The  experimental  M  =  0.6  contour  is  a  nearly  uniform,  radial 
one  which  explains  the  agreement  between  the  one-dimensional  analysis  and 
the  wall  measurements  of  Fig.  iV.11  in  the  upstream  region. 


In  view  of  the 
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transonic  approximati ons  that  have  been  made,  the  series  solution  developed 
in  the  preceding  chapter  would  be  expected  to  accurately  predict  the  throat 
flowfields  only  for  regions  near  the  sonic  line.  Thus,  the  agreement 
between  this  analysis  and  the  experimental  data  for  Mach  numbers  up  to  1.4 
through  such  a  wide  region  of  the  throat  is  somewhat  surprising,  particularly 
for  the  small  wall  radius  of  curvature  of  this  model. 

In  Fig.  IV. 13  the  expansion  solution  results  are  compared  to  the  splitter 
plate  flowfield  measurements  for  an  annular  nozzle  whose  inner  wall  radius  of 
curvature  is  twice  that  of  the  outer  wall.  This  is  the  configuration  drawn 
in  Fig.  IV. 10  for  which  the  centers  of  curvature  of  both  wall  contours  and, 
therefore,  the  throat  all  lie  along  the  Z  =  0  plane  so  that  the  one-dimensional 
flow  direction  is  parallel  to  the  axis  of  symmetry,  6=0.  Since  the  maximum 
radius  of  the  centerbody  is  0.375  units  in  these  R-Z  coordinates,  the  throat 
half-height  is  reduced  to  d  =  0.625  resulting  in  dimensionless  radii  of 
curvature  of  R  =  3.2  and  R  =1.6  for  the  inner  and  outer  walls, 
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respectively.  The  dimensionless  distance  to  the  inner  throat  wall  location 
is  y.  =  0.6.  The  data  points  have  been  obtained  in  a  manner  similar  to  that 
described  for  the  conventional  axisymmetric  geometry  of  Fig.  IV. 12  except 
that  for  this  and  the  other  two  annular  configurations  investigated,  there 
are  five  rows  of  pressure  taps  at  a  constant  radial  coordinate  in  addition 
to  the  inner  and  outer  walls:  R  =  7/16,  9/16,  11/16,  13/16,  and  15/16.  It 
is  seen  that  the  theoretical  results  correspond  very  closely  to  the  measure¬ 
ments  except,  perhaps,  near  the  tangency  point  between  the  conically  conver¬ 
gent  section  and  the  circular  arc  contour  along  the  outer  wall,  Z  =  -  0.6. 
Otherwise,  the  agreement  is  very  good  through  the  entire  throat  region  from 
M  =  0.6  to  M  =  1.4. 


Figure  IV. 14  presents  the  comparison  between  the  theoretical  and  experi¬ 
mental  results  for  the  annular  configuration  in  which  the  centerbody  has  been 
translated  downstream  relative  to  the  outer  wall  contour  by  a  distance  of 
12.7  mm  or  0.5  units  in  the  R-Z  coordinate  system  of  the  figure.  For  this 
geometry  the  main  flow  direction  is  inclined  at  an  angle  of  approximately 
5.45°  (0.095  radians)  away  from  the  axis  of  symmetry.  Also,  the  throat  half¬ 
height  is  larger  than  for  the  preceding  case  in  which  the  centers  of  curva¬ 
ture  of  both  walls  were  on  the  Z  =  0  plane  so  that  the  radii  of  curvature 
and  distance  to  the  inner  throat  wall  location  in  the  local  x-y  coordinates 
are  all  reduced:  R  =  3.01,  R  =  1.50,  and  y  =  0.51.  Since  this  is  the 
first  example  for  which  the  location  of  the  minimum  area  cross  section  is  not 
obvious,  it  has  been  plotted  as  the  dashed  line  in  the  figure.  Also  shown 
by  the  dotted  line  is  the  cross  section  of  minimum  distance  between  the  two 
contours  which  joins  their  centers  of  curvature.  As  previously  mentioned,  the 
minimum  area  throat  cross  section  for  inclined  nozzles  is  nearer  the  axis  of 
symmetry  than  is  the  cross  section  of  minimum  distance,  and,  as  can  be  seen  in 
the  figure,  the  two  are  not  necessarily  parallel.  Even  for  the  relatively 
small  inclination  angle  of  this  configuration,  there  is  a  significant  distance 
separating  the  two  cross  sections  because  of  the  relative  proximity  of  the 
throat  to  the  axis  of  symmetry.  In  comparing  the  analytical  and  experimental 
results  in  this  figure,  it  is  noted  that  the  correspondence  is  quite  close 
except  in  the  upstream  regions  near  the  two  walls  where  the  circular  arc  con¬ 
tours  of  the  experimental  models  end.  Along  the  centerbody,  the  circular 
inner  contour  meets  the  upstream  cylindrical  section  at  Z  =  -  0.2,  while  the 
circular  arc  outer  contour  is  tangent  to  the  conical  inlet  at  Z  =  -  0.6. 


Elsewhere  in  the  flowfield,  the  agreement  between  the  theoretical  and  experi 
mental  iso-Mach  curves  from  M  =  0.6  to  M  =  1.4  is  quite  good. 

The  final  experimental  geometry  investigated  is  shown  in  Fig.  IV. 15.  I 
this  case  the  centerbody  was  displaced  upstream  relative  to  the  outer  nozzle 
wall  by  12.7  mm  or  0.5  units  in  the  R-Z  coordinate  system.  The  one¬ 
dimensional  flow  direction  is  inclined  toward  the  axis  of  symmetry  at  an 
angle  of  5.45°,  and  the  cross  sections  of  minimum  area  and  minimum  distance 
between  the  contours  are  identical  to  those  of  the  preceding  configuration 
as  reflected  through  the  Z  =  0  plane.  The  radii  of  curvature  and  the  dis¬ 
tance  to  the  inner  nozzle  wall  also  have  the  same  values  as  the  last  example 
R  =  3.01,  R  =  1.50,  and  y  =  0.51.  For  this  case  the  theoretical  and 
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measured  Mach  number  contours  are  compared  only  for  Mach  numbers  up  to 
M  =  1.2.  For  higher  Mach  numbers,  the  measurements  were  disturbed  by  the 
presence  of  a  detached,  oblique  shock  wave  emanating  from  the  inner  wall 
near  Z  =  0.2  where  the  supersonic  flow  along  the  centerbody  is  turned  paral¬ 
lel  to  the  downstream  cylindrical  section.  The  series  solution  results 
again  agree  quite  well  with  the  experimental  data  except  in  the  upstream 
region  along  the  outer  wall  near  the  Z  =  -  0.6  tangency  point. 

Cross-plots  comparing  the  experimental  wall  Mach  numbers  with  those 
determined  from  the  expansion  solution  and  from  isentropic,  one-dimensional 
analysis  could  also  be  presented  for  the  three  annular  configurations  as  was 
done  for  the  conventional,  axisymmetric  nozzle.  Fig.  IV. 11.  However,  as  can 
be  seen  by  studying  Figs.  IV.13-IV.15,  the  resulting  figures  and  conclusions 
would  all  be  similar  to  those  for  the  Laval  nozzle.  Near  the  throat  and 
downstream  of  it,  the  series  solution  agrees  well  with  the  measurements  of 
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the  Mach  number  along  the  walls,  while  the  one-dimensional  results  are  in 
error.  However,  in  the  upstream  regions  where  the  circular  arc  wall  con¬ 
tours  of  the  models  make  transitions  to  either  conical  or  cylindrical 
sections,  the  agreement  between  the  expansion  solution  and  the  data  is 
understandably  not  as  good. 

Without  exception,  it  is  noted  that  in  regions  of  the  flowfields 
unaffected  by  the  wall  contour  transitions,  the  shapes  of  the  experimental 
and  theoretical  Mach  number  contours  are  similar  but  the  measurements  lie 
slightly  downstream  of  the  analytical  results.  This  observation  applies  not 
only  to  the  comparisons  made  in  this  section.  Figs.  IV .  1 2- IV . 1 5,  but  also 
to  the  comparisons  made  with  existing  data  sets  in  Chapter  III,  Figs.  III. 10 
and  III. 11. 

One  cause  of  this  behavior  results  from  the  fact  that  the  analysis  is 
an  inviscid  one  while  the  effect  of  wall  friction  causes  a  shift  in  the 
actual  sonic  line  downstream  from  the  corresponding  inviscid  one.  This  can 
be  shown  by  means  of  the  following  argument.  For  a  perfect  gas  with  constant 
specific  heats,  one  of  the  differential  equations  governing  generalized  one¬ 
dimensional  flow  with  friction  and  area  change  may  be  written  as  [90], 

(i-m2)  ir  =  %-  r  >  (IV_5) 

where  u  is  the  one-dimensional  velocity,  dF  is  the  elemental  wall  shear  force 
opposing  the  flow,  p  is  the  static  pressure,  and  A  is  the  flow  cross-sectional 
area.  At  the  sonic  line  the  left  hand  side  of  this  equation  vanishes,  resulting 
in, 
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fn,r  me  wail  snear  forte  and  the  static  pressure  are  inherently  positive,  dA 
"jst  also  re  a  positive  quantity.  It  is  therefore  concluded  that  the  effect 
o*  friction  is  to  shift  the  sonic  line  downstream  into  the  diverging  part  of 
the  supersonic  nozzle  for  which  dA/dx  0. 

A  secono  reason  that  the  experimentally  determined  iso-Mach  curves  are 
generally  slightly  downstream  from  the  theoretical  predictions  is  related  to 
measurement  errors  introduced  by  the  static  holes.  Since  the  pressure  cali¬ 
brations  were  performed  in  a  no-flow  test  in  which  the  transducers  were 
directly  connected  to  a  pressure  source,  the  accuracy  estimate  made  in  the 
previous  section  considered  only  the  errors  resulting  from  the  transducers, 
signal  conditioning  equipment,  calibration  procedure  and  equipment,  etc. 
However,  with  flow  across  the  holes  an  additional  error  is  present  since  it 
is  well  known  [91-94]  tnat  the  pressure  measured  by  a  static  hole  is  not 
necessarily  identical  to  the  static  pressure  in  the  flowfield  at  that  point. 

The  pressure  tap  configuration  used  in  this  study  is  shown  in  Fig.  IV. 16. 
It  is  very  similar  to  the  one  used  by  Shaw  [92]  in  his  study  of  pressure  tap 
errors  with  the  inside  diameter  of  the  tap  twice  that  of  the  hole.  Also,  in 
the  present  experiments  the  hole  depth  to  diameter  ratio,  i/d.  ,  was  nominally 
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equal  to  two.  Using  dimensional  analysis  and  the  results  of  extensive  exper¬ 
iments  for  incompressible  turbulent  flow,  Shaw  correlated  the  static  pressure 
errors  as  a  function  of  the  hole  depth-to-diameter  ratio  and  the  Reynolds 
number  based  on  the  hole  diameter  and  friction  velocity, 

f2  =  f[Re+,Z/dhj  ( I V- 7 ) 


where 


i  iv-b) 


f  e 

is  the  wall  shear  stress  and  .  is  the  kinematic  viscosity  or  the  fluid. 
Shaw  found  that  for  deep  holes,  .'/d  1.5,  the  pressure  error  is  positive 

and  is  independent  of  the  hole  depth,  '.  The  dimensionless  error,  :.p/ :  , 
increases  from  zero  at  Re+  =  0  until  it  reaches  an  asymptotic  value  of  2.75 
at  Re*  --  800,  above  which  the  error  is  constant. 

Using  a  somewhat  different  pressure  tap  configuration,  Livesey,  et  al. 
[93]  performed  a  similar  set  of  experiments  and  analysis  of  static  hole 
errors.  Tne  results  were  correlated  in  the  same  form  as  tiiose  of  Shaw,  but 
for  deep  holes  it  was  found  that  the  positive  error  does  not  reach  an 
asymptotic  limit  out  instead  increases  monotonical ly  and  nearly  linearly  for 
all  hole  Reynolds  numbers.  Below  Re+  =  1000  the  error  predicted  by  the 
correlation  of  Livesey,  et  al.  is  lower  than  that  found  by  Shaw.  At 
Re+  =  1000  the  correlations  intersect  so  that  above  this  point  the  monotonic- 
ally  increasing  error  predicted  by  Livesey,  et  al .  exceeds  the  limiting  value 
of  Shaw. 

It  is  important  to  note  that  the  positive  static  pressure  errors  found 
by  both  of  these  investigators  are  consistent  with  the  trends  exhibited  in 
Figs.  III. 10,  III. 11,  and  IV.12-IV.15.  For  these  cases  the  Mach  number 
measured  at  a  given  point  is  lower  than  the  corresponding  predicted  value 
so  that  the  static  pressure  is  higher.  Therefore,  taking  into  account  the 
predicted  static  hole  errors  puts  the  data  and  theory  in  closer  agreement. 

By  using  the  boundary  layer  analysis  developed  in  [88]  to  predict  the  wall 
shear  stress,  it  is  found  that  the  hole  Reynolds  numbers  for  the  experiments 
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of  this  investigation  ranged  from  Re+  =  870  at  M  =  0.6  Lu  Rc+  =  1120  at 
M  =  1.4.  For  these  values,  the  errors  predicted  by  Shaw  and  Livesey,  et  al. 
agree  quite  closely  so  that  Shaw's  limit  of  .Lp/r0  -  2.75  was  used.  For 
M  =  0.6  this  corresponds  to  a  static  pressure  error  of  approximately  +  0.2' 
while  for  K  =  1.0  and  M  =  1.4  the  errors  are  calculated  as  +  0.5%  and  +  O.9., 
respectively.  The  percent  root-mean- square  difference  between  the  theoreti¬ 
cal  and  measured  values  of  the  static  pressure  for  the  data  presented  in 
Figs.  IV. 12- IV. 15  is  somewhat  higher  than  these  predictions  for  the  static 
pressure  error.  For  the  axisymmetric,  Laval  nozzle  of  Fig.  IV. 12  the  per¬ 
centage  difference  is  found  to  be  (Ap/p.)  =  3.5%,  while  for  the  annular 

configurations  of  Figs.  IV. 13,  IV. 14,  and  IV. 15  the  differences  are  2.5%, 
2.6%,  and  3.3%,  respectively.  It  should  be  remembered,  however,  that  the 
investigations  of  both  Shaw  and  Livesey,  et  al.  considered  only  incompress¬ 
ible  flow.  Rayle  [91],  as  reported  by  Chue  [94],  found  that  for  transonic 
speeds  the  static  pressure  error  is  two  to  three  times  that  for  incompress¬ 
ible  flow  for  the  O.ol  mm  diameter  static  holes  used  in  this  investigation. 
Thus,  the  static  hole  errors  may  be  as  large  as  2-3%  for  these  experiments. 

It  is  felt  that  the  major  reasons  for  the  downstream  shift  of  the 
experimental  constant  Mach  number  curves  from  the  corresponding  analytical 
ones  are  the  two  just  discussed:  wall  friction,  which  is  not  accounted  for 
in  the  theoretical  model,  and  errors  in  the  static  pressure  measurements. 
Other  effects  such  as  heat  transfer,  condensation,  particulates  suspended 
in  the  flow,  etc.  are  thought  to  be  negligible.  It  should  be  emphasized 
that  the  discrepancies  are  small  and,  for  the  most  part,  the  agreement 
between  the  measurements  and  the  theory  is  quite  good  through  a  surprisingly 
wide  region  of  the  throat  flowfields. 


Figure  IV. 1  Photograph  of  full-section,  conventional 
axisymmetric  nozzle 
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Figure  IV. 4  View  downstream  through  test  section  of  assembled 
half-section,  annular  nozzle 
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Figure  IV. 5  View  upstream  through  measurement  region  of  assembled, 
half-section,  annular  nozzle 


Figure  IV. 6  Partially  assembled,  annular  nozzle  test  apparatus 


Figure  IV. 7  View  of  annular  nozzle  through  side  window  of 
viewing  chamber 


Figure  IV. 8  Overall  view  of  fully  assembled,  annular  nozzle 
test  apparatus 
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Figure  IV. 10  Half-section  drawing  of  annular  axisynmetric  experimental 
configuration;  all  dimensions  are  in  millimeters 
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Figure  IV. 15  Comparison  of  constant  Mach  number  contours  from  series 
expansion  solution  with  experimental  data  for  annular 
nozzle  with  centerbody  center  of  ^urvature  along 
Z  =  -  0.5  plane;  Re2d  =  2.08  x  106  for  experiments 
(Flagged  data  points  obtained  by  extrapolation] 
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Flow 


Figure  IV. 16  Half-section  drawing  of  static  hole-pressure  tap 
arrangement;  all  dimensions  are  in  millimeters 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 

As  a  result  of  the  integrated  theoretical  and  experimental  investigation 
of  nozzle  throat  flowfields  which  has  been  reported  herein,  the  following 
conclusions  may  be  drawn: 

(1)  The  approximate  series  expansion  solution  which  has  been  developed 
provides  a  direct  means  of  analyzing  the  transonic  flow  in  the  throat  region 
of  a  wide  variety  of  two-dimensional  nozzle  configurations  including  axi- 
symmetric,  annular,  and  planar  nozzles.  For  annular  nozzles,  the  throat  may 
be  inclined  with  respect  to  the  axis  of  symmetry,  and  for  both  planar  and 
annular  geometries  the  two  bounding  walls  may  have  unequal  radii  of  curva¬ 
ture  at  the  throat. 

(2)  Based  on  the  series  of  parametric  studies  which  was  performed,  it 
is  suggested  that  third  order  solutions  be  employed  together  with  n  =  2, 
where  n  is  the  parameter  in  the  expansion  variable  definition,  Eq.  (III-22). 

(3)  From  the  parametric  studies  and  from  comparisons  between  the  series 
solution  results  and  experimental  data,  it  is  concluded  that  the  expansion 
solution  accurately  describes  the  transonic  flowfield  through  a  significant 
region  of  the  throat  for  nozzles  whose  wall  radius  of  curvature  is  of  the 
order  of  the  throat  half-height  or  larger.  In  addition,  for  an  axisymmetric 
nozzle  with  the  small  dimensionless  wall  radius  of  curvature  of  R  =  0.625, 

C 

it  was  found  that  the  series  solution  is  accurate  from  the  sonic  line  up  to 
a  Mach  number  of  approximately  1.4. 

(4)  An  accurate,  slightly  supersonic  starting  line  for  initiating  the 
hyperbolic  computations  for  the  supersonic  portion  of  the  flowfield  is  the 
constant  Mach  number  contour  emanating  from  one  of  the  throat  wall  locations. 
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This  initial  value  line  is  recommended  because  the  wall  boundary  conditions 
are  exactly  satisfied  only  at  the  throat.  It  may  be  more  convenient, 
however,  to  use  other  starting  lines  for  compatibility  with  the  scheme  used 
to  analyze  the  supersonic  flowfield. 

(5)  The  numerical  implementation  of  the  expansion  solution  is 
extremely  fast  and  reliable  since  it  simply  consists  of  the  straightforward 
evaluation  of  a  series  of  constants  and  functions.  Typical  problems  require 
less  than  two  seconds  of  central  processor  (CPU)  time  on  the  University  of 
Illinois  CDC  Cyber  175  digital  computer.  This  feature  makes  feasible  para¬ 
metric  studies  and  iterative  calculations  as  might  be  necessitated  in  the 
nozzle  design  situation. 

(6)  For  applications  in  which  a  very  high  degree  of  accuracy  is 
required,  it  may  be  advantageous  to  employ  the  approximate  series  solution 
as  the  initial  condition  for  one  of  the  time-dependent  numerical  techniques. 
This  scheme  may  speed  the  convergence  of  the  time-dependent  methods  to  the 
steady  state  asymptotic  limit. 

(7)  Based  on  the  comparison  of  the  wall  pressure  measurements  from  the 
full-  and  half-section  models  of  the  conventional,  axisymmetric  nozzle,  it 
is  concluded  that  the  splitter  plate  technique  is  valid  for  obtaining 
flowfield  pressure  distributions  for  these  high  Reynolds  number,  accelerated 
flows.  The  boundary  layer  on  the  splitter  plate  is  extremely  thin  so  that 
it  has  little  effect  on  the  flowfield  pressure  measurements. 

(8)  For  the  axisymmetric,  Laval  nozzle  and  the  three  annular  configura¬ 
tions  investigated  here,  the  agreement  between  the  theoretical  and  experimental 
results  for  the  constant  Mach  number  contours  is  quite  good  through  a 
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significant  region  of  the  throat.  Only  near  the  points  where  the  circular 
arc  contours  of  the  experimental  models  make  transitions  to  either  the 
conical  inlet  or  cylindrical  approach  sections  is  there  a  noticeable  dis¬ 
crepancy  between  the  analysis  and  the  measurements. 

(9)  The  small  downstream  shift  of  the  experimental  iso-Mach  curves 
from  the  corresponding  theoretical  ones  can  be  attributed  to  wall  friction, 
which  is  neglected  in  the  theory,  and  to  pressure  measurement  errors 
related  to  the  static  pressure  taps. 

The  understanding  of  the  topic  of  internal  transonic  flowfields  is  by 
no  means  complete.  To  further  enhance  the  state  of  knowledge  in  this  area, 
the  following  studies  are  recommended: 

(1)  Through  the  use  of  alternate  expansion  techniques  it  may  be  possi¬ 
ble  to  obtain  solutions  which  are  valid  for  nozzles  whose  throat  wall  radius 
of  curvature  is  small  compared  to  the  throat  >->*lf- height.  Alternate  expan¬ 
sions  should  also  be  investigated  for  the  case  of  highly  inclined,  annular 
nozzles,  particularly  for  the  radial  flow  case  which  has  practical  applications. 

(2)  The  possibility  of  extending  the  series  solution  technique  to  uncon¬ 
ventional,  three-dimensional,  supersonic  nozzles  should  also  be  investigated. 

It  is  felt,  however,  that  the  most  likely  candidates  for  analyzing  the  flow- 
fields  in  this  class  of  nozzles  are  numerical  methods. 

(3)  In  order  to  investigate  the  effect  of  the  Reynolds  number  on  these 
flows,  experimentation  over  a  wider  range  of  pressures  is  needed.  Additional 
experiments  are  also  required  for  nozzles  having  a  small  throat  wall  radius 
of  curvature  and  for  highly  inclined  annular  nozzles. 

(4)  For  further  flowfield  information,  such  as  detailed  measurements  of 
the  velocity  components,  a  laser  doppler  velocimeter  system  may  be  employed. 
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The  theoretical  solution  and  experimental  static  pressure  measurements  of 
the  present  investigation  provide  a  basis  for  setting  up  and  checking  out 
such  a  system. 

(5)  An  area  which  seemingly  has  received  very  little  attention  is  that 
of  converging-diverging  nozzles  operating  in  the  unchoked,  Taylor  regime 
for  which  the  operating  pressure  ratio  is  important.  An  important  applica¬ 
tion  of  this  problem  is  jet  engine  inlets  for  subsonic  aircraft. 

(6)  Ir,  a  related  area,  it  appears  that  additional  work  is  needed  to 
better  understand  static  pressure  tap  errors  which  are  incurred  in  compress¬ 
ible  flow  measurements. 
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APPENDIX  A.  INVESTIGATION  OF  EQUATIONS  FOR 
CIRCULAR  ARC  CONTOURS 


Assuming  that  the  bounding  walls  of  the  annular  nozzle  are  circular- 
arcs,  their  equations  in  the  x-y  coordinate  system  may  be  written  as 


(x-a  )2  +  (y-b  )2  =  R  2 

1  1  Cl 

(A- 1  ) 

and  (x-a  )2  +  (y-b  )2  =  R  2  , 

o  o  CO 

(A-2) 

where  (a  ,b.  )  and  (a  ,b  )  are  the  coordinates  of  the  centers  and  R  and  R 

11  o  o  c  1  c  o 

are  the  radii.  Fig.  A.l.  In  the  Maclaurin  series  expansions  of  these  con¬ 
tour  equations, 


y  =  g(x)=y;  +  g 1 (0 ) x  +  g" ( 0 )  +  g'"  (0)  ~  + 
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y  =  h  ( x )  =  yQ  +  h '  (0  )x  +  h"(0)  ~r  +  h"1  (0)  +  ... 

the  coefficients  can  be  calculated  as 
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From  the  definition  of  the  expansion  parameter,  Eq.  (Ill - 20 ) ,  it  is  known 
that  the  radii  of  curvature,  R  and  R  ,are  (c-1)-  However,  as  can  be 

c  i  co 

seen  from  Eqs.  (A-5)  and  (A-6),  estimates  for  the  x  coordinates  of  the 
centers,  a.  and  a  ,  must  also  be  obtained  in  order  to  complete  the  order  of 
magnitude  analysis  of  the  Maclaurin  series  coefficients. 

To  do  this,  the  cross-sectional  flow  area  in  the  nozzle  throat  in  x-y 
coordinates  is  first  written  as 

rh  < x  >  r  ch  <  *  > 

A(x)  =  2tt  ~  dy  =  2tt  (x  sind+y  cos6)dy  ,  (A-7) 

x)  t  x) 

where  the  transformations  in  Eqs.  ( 1 1 1- 3)  and  ( 1 1 1 - 4 )  have  been  used.  Per¬ 
forming  the  integration  and  using  the  expansions  in  Eqs.  (A-3)  and  (A-4), 
the  first  few  terms  in  the  expression  for  the  area  are 

A|^)_  =  cosd  j/_y2j  +  sin3(y  -y  )  +  cos6(y  h1 (0)-y.  g'(0))lx 

CTT  C  ,  °  1  j  L  °  1  °  1 

( A- 8 ) 

+  sin8(h'(0)-g'(0))  +  [yoh"(0)-y.  g"(0)  +  h,(0)2-g'(0)2]  x2 

+  . . . 

Since  by  definition  the  area  is  minimum  at  the  throat  plane,  x  =  0,  dA/dx 
must  vanish  there.  This  leads  to  the  relation 

tanB  =  y.  g'(0)  -  yoh'(0)  ,  cosB  t  0r  (A-9) 

since  yo~yi  -  1.  Defining  an  average,  symmetric  throat  slope,  j'(0),  and  an 
average  y  location  in  the  nozzle  throat,  y,  by 

2yj'(0)  =  y.  g* (0)  -  yo h ■ (0)  ,  (A-10) 

fThe  radial  flow  case,  $  =  ir/2,  was  previously  ruled  out,  Eq.  ( 1 1 1-28 ) . 
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tq.  (A-9)  yields  the  following  relation  for  j'(0), 

j*  (0)  *  \  .  (A- 1 1 ) 

In  view  of  the  previous  restriction  placed  on  tanE/y,  Eq.  (Ill- 30 ) ,  the 
throat  slope  quantities,  j'(0),  g‘(0),  and  h'(0),  may  be  estimated  as 
(>.  ).  Combining  this  information  with  the  expressions  for  g *  ( 0 )  and 

h ' (0 )  from  tqs.  (A- 5 )  and  (A-6)  and  the  fact  that  the  radii  R  and  R  are 

c  i  co 

- (c  1 )  results  in  the  conclusion  that  the  x  coordinates  of  the  centers  of 

l/2 

curvature,  a,  and  ao ,  are  '( 0  )  under  the  present  assumptions.  The  order 

of  magnitude  estimates  for  the  Maclaurin  series  coefficients  for  circular 
arc  boundaries  can  then  be  completed  as, 

g'(0),  h* (0)  =  0(e3/2)  9" (0) ,  h"(0)  =  o(±) 

,  (A- 12) 

g"'  (0),  h"1  (o)  =  o{z  2 )  glv  (0),  h,v  (0)  =  o(c3) 

•  • 

•  • 

•  • 

Although  the  details  of  the  boundary  condition  evaluations  have  been 
presented  here  only  for  circular  contours,  the  results  are  also  valid  for 
other  contours  whose  Maclaurin  series  coefficients  have  similar  orders  of 
magnitude.  Thus,  the  solutions  developed  in  Chapter  III  are  expected  to 
apply  to  other  conic  section  contours  such  as  parabolic,  hyperbolic,  and 
elliptic  arcs.  It  should  also  be  noted  that  the  boundary  conditions  at 
straight  walls  are  identically  satisfied  since  the  wall  slope  is  exactly 
matched  and  all  higher  derivatives  vanish  for  this  case. 


gure  A.1  Configuration  and  nomenclature  for 
circular  arc  contour  equations 
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APPENDIX  B.  SOLUTION  SUMMARY 

Referring  to  Fig.  III.l,  the  only  quantities  which  need  to  be  specified 
in  order  to  initiate  the  solution  are:  y*  the  ratio  of  specific  heats;  r, 
the  parameter  in  the  expansion  variable,  Eq.  (III-22);  N,  the  number  of 
terms  of  the  series  to  be  included;  and  the  equations  for  the  inner  and 
outer  wall  contours,  G(R,Z)=0  and  H(R,Z)=0.  With  the  geometrical  configura¬ 
tion  given,  the  following  quantities  may  be  determined:  the  minimum  area 
cross-section,  i.e.,  the  throat  wall  points  (Z.  ,R.  )  and  {Zq,Ro);  d,  the 
throat  half- height;  Z*,  the  Z-location  of  the  x-y  origin;  and  6,  the  angle 
of  inclination  of  the  x-axis  with  respect  to  the  symmetry  axis.  Utilizing 
the  coordinate  transformations  listed  in  Eqs.  ( 1 1 1 - 3 )  and  (III-4),  the 
parameters  y(  ,  y  ,  g ' (0 ) ,  h'(0),  g"(0),  and  h"(0)  in  the  x-y  coordinate 
system  can  be  evaluated;  y  =  g(x)  and  y  =  h(x)  are  the  transformed  equations 
of  the  boundaries. 

In  order  to  carry  out  the  initialization  of  the  solution  as  well  as 
its  further  evaluation,  a  library  of  FORTRAN  routines  has  been  developed 
and  used.  The  throat  plane  is  located  numerically  using  a  pivoting  scheme 
whereby  the  location  of  the  inner  or  outer  wall  point  is  alternately  fixed 
and  a  sweep  is  made  across  the  opposite  contour  until  the  minimum  area  sec¬ 
tion  for  the  given  fixed  point  is  found.  This  process  is  continued  until 
the  fractional  change  in  area  between  iterations  is  less  than  a  prespecified 
convergence  value,  e.g.,  10’10.  This  method  has  been  thoroughly  checked  by 
substituting  results  obtained  with  it  into  the  Lagrange  multiplier  formula¬ 
tion  of  the  problem  of  minimizing  the  flow  cross-sectional  area  in  an 
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annular  nozzle  with  circular  arc  walls.  Since  the  resulting  equations  have 
been  found  to  be  satisfied  to  a  high  degree  of  numerical  accuracy,  e.g., 

__  5 

residuals  on  the  order  of  10  or  smaller,  it  is  concluded  that  this  method 
can  be  used  to  locate  accurately  and  reliably  the  throat  plane.  The  contour 
derivatives  g’(0),  h'(0),  g"(0),  and  h"(0)  are  also  evaluated  numerically 
using  second  order  accurate,  centered  differences  and  a  nodal  spacing  of 

—  3  -6 

Ax  =  10  .  The  truncation  error  is  therefore  on  the  order  of  10  which  is 

sufficiently  accurate  for  the  present  purposes. 

Once  the  quantities  mentioned  above  have  been  found,  all  of  the  flow- 
field  properties  of  interest  can  be  determined  by  the  straightforward 
evaluation  of  a  series  of  constants  and  functions.  Given  the  point  (Z,R) 
at  which  the  information  is  desired,  transformation  Eqs.  ( 1 1 1- 3 )  and  (III-4) 
are  first  employed  to  find  the  corresponding  (x,y)  coordinates  and  the  fol¬ 
lowing  parameters  are  then  evaluated, 


e  = 
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The  transonic  perturbation  velocity  components 


(u3,v3)  are  then  found  as  follows. 


(B-l ) 
(B-2) 
(B-3) 
(B-4) 
(B-5 ) 
(B-6) 

(Uj  ,Vj ),  (u2 ,v2 ) ,  and 
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First  Order  Solution 


The  u  and  v]  velocity  components  from  the  first  order  solution  are 
given  by, 

Uj  (z,y)  =  A]  (y)  +  B0  +  8^  (B-7) 

Vj (z,y)  =  Ag(y)  +  Aj (y)z  ,  (B-8) 

where  the  constants  are  determined  in  the  following  order. 
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and  the  functions  of  y  are  determined  as, 

A;(y)  =  +J- 

Aj  (y)  =  }  Bjy2  +  B2lny 
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The  discharge  coefficient  constant  C  is  found  as. 


CD4  ■  J  B0B!  -  16  B>> 


D6 


-  1  R2  +  1  R2 

2  o 

4  2 

=  B  B  - 

4  b2 

0  2 

2  2 

1  n4  i 

t 

6  6 

24  Bi 

y  -y. 

J  o  •'i 

+  c. 


D4 


4  4 

y  -y. 

J  o  Ji 


L  r  22 
+  c  y  -y 

D5  Jo 


+  1  R2 
+  2  B2 


y„  i  ^y. 


(B-19) 
(B-20) 
( B-21 ) 


2-y?U«y, 


+  |  b2b 

4  1  2 


y4Zny  -y4  Zny. 

O  O  I  J 


+  c, 


D6 


y2olnyQ-yhny. 


(B-22) 


Second  Order  Solution 

The  u2  and  v2  velocity  components  from  the  second  order  solution  are 
given  by 


ua(z,y)  =  [Cj  (y)  +  D0]  +  (2C2(y)  +  DJz  +  D2z2  (B-23) 

v2(z,y)  =  Cg(y)  +  C;(y)z  +  C^{y)z2  ,  (B-24) 


where  the  constants  are  determined  in  the  order, 
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C  =  C  +  C  +  C 

D2  D7  D8  D9 


( B-89 ) 


Third  Order  Solution 

The  u3  and  v3  velocity  components  from  the  third  order  solution  are 
given  by, 

u3(z,y)  =  [E!  (y)+F0j  +  [2E2  (y)+Fxj  2  +  [3E3(y)+F2lz2  +  F3z3  (B-90) 

v3(z,y)  =  E'(y)  +  E;(y)z  +  E^(y)z2  +  E'(y)z3  ,  (B-91) 

where  the  constants  are  determined  in  the  following  order, 

F.  ■  hi  [j6  °<A  +  +  l  i„»!!"y, 

*  “i.^y.  ♦  “„]  <B-92> 

,  1 

+  h. 


Fs  =  h2n  +  h2n 


\  +  B2Z^o  +  B0 


-  BlVo  +  K[ln*o)2  +  1D14^^0  +  D16Z^o  +  °0 


T6  D13y!  +  D20^ 


+  h. 


Dlz*o  +  2Dn^o  +  D1 


rti 

2  : 


|bV 

4  i^o 


+  B6  +  W”**  +  B!B2  - 


(B-93) 


F6  =  h2nBj  +  h1D2  +  h2^D12y2  +  2Du^yo  +  D, 


2  B2‘ 

B  -  — 
1  2 

y„ 


M 

g 

1  ♦  b,  *  *  B,v  4 

- 1 

>>° 

_ 1 

f,  =  M 

7  2  2 


fel1 

h 

r  b  i 

B2-  ~ 

l  4 

2 

1  y2 

Jo 

(B-94) 


(B-95) 


143 


F  =  g  4-^  D  y  +  D  y  -  B  y  +  B'  1  '•  y  '*  +  |  D  y  •  y 
6  yl  16  13Ji  20Ji  1  1  y\  2  ,  y,  |  4  !  4~i  y, 


+  D  >;y  +  D  ■ 

16  0 


{ B- ?6 ) 


F  -  g  •  +  g  g  i  i  B  y  +  B  •  y  +  B  +  g  1 D  y  +  20  ■.  y  +  0 

9  d2  y2  ;  2  2  0  yl  12J.  11  yi  1 


+  n  il-  n  U4  +  n  V2-  •  R  W  +  R2  I  ‘  »  V  )2  ■  1 


2  . 


g  i  -r>-  D  y  +  D  y  -  •  S  y  +  B  !  .  >.  y  +  -,  D  v  -y 
y2  16  13Ji  2CJi  1  l^i  2  i  yi  i  4  M-'i  g 

+  D  3>:y  +  D 

i  a  J  i  n 


f  \  + 1  1  1  3  3  2 

-  -  -  -  -V-i  9  !3t  B  y  +  B  +  B  B  \y 

It.  i  Oj  ^  2  !  yl  !  4  1  yi  6  1  2 


VI 


+  B,B2  -  T' 

y. 


(B-97) 


i  ^,+r- 


F„=gnB  +  q  D  +  q  D  y  +2D  b:  y  +  D  -  ■4-=—,  q 
10  1  yl  2  y2  |_  12Ji  11  ij  ^  2  j  y 


111 

4 


I  ntf  *  8.  *  B,  B2  Z„y,  ♦  6,B2  -  4 

y, 


V 


y, 


( B-98) 


F  =  g  D  - 
11  y2  2  4 


2  8  2 

B  -  — 
1  2 

y, 


F12  =  2(y-l)B1D12  +  2(y-1)b'D2  +  (y-l)B*  +  4D*  +  4bJD2 


F  ■ 


1  3 


F3  = 


2  2 

y  -y 

J  o  ■'i 

c  1  c  2  c 

F7V  2  F12^o-F 


1  3 


4Bl*o 


( B-99 ) 

(B-100) 

(B-101) 

(B-102) 


F14  =iF12  +4B1F3 


(B-103) 


144 


F15  =  3Bj  F3  +  (Y+2)b|d2  +  9BlFM  +  6D2D12  +  3(Y+2)B>12  +  (2^ 


+  (y-DB.Dj. 


(B-104) 


F,»  *  f1!1]3,8!  -  6i°2 
F,,  =  l5?1)8,8'  +  682D1  1 
Fi«  '  (1rlBiB2 


(B-105) 


(B-106) 


(B-107) 


F19  »  6B2F3  +  2(Y+2)BlB2D2  +  18BiF13  ♦  1 20a  D, ,  +  2(2Y+1)B>ii 

+  2(Y-1)B2D12  +  (Y-1)BiD14  +  3{Y-1)B^B2  (B-108) 


F2o  =  P^KB2  +  6BX  +  6B2°i2  +  2(Y-1)B0D12  -  2(y-DBi  Dx 

+  M8.0!.  *  2(v-')BsDi  +  (y-I)bX  ♦  (v-1)B,B2D2 


+  (2y+1)BJdi  +  2(Y-1)B0b'  +  6BqF3  +  60^  +  ^0^D2  (B-109) 


1  i  ^ny  F  r,  , 

f22  =  i  Fi5y0  +  p16  +  p17  -y- '  —  +  F* 9 1  yj*y0  -  i  yc 

y  L 


+  2  F20yo 


(B-no) 


c  1  c  3  lnyi  p18  r  1  7  1 

F23  =  4  F15yi  +  F16  +  F17  "y - —  +  F19  2  U*i  '  4  yi 

Ji  V .  L  J 


+  2  F20yi 


F10yo  ~F6yi  ’F23yo+F22yi  |yi  yo 


2  2 
y  -y 

Jo 


(B-ni) 


(B-112) 


145 


F  y  -F  -F  y 

bJo  21  22}o 


r 

2 

(B-113) 

F2  4 

s  7  F„  - 

l  r„  *  3B, 

F 

2 

(B-114) 

F2  5 

’  ?  F»  - 

5  F„ 

(B-115) 

F26 

■  3<F„  * 

21  B,F,S  * 

20  2  +  -!-  0  D 

1  2  4  2  1  3 

I5_1+4')R2 
*  l  “  )Bi 

°,3  * 

*  \  BX 

l  2 

+  (y-1)B1D2i 

(B-116) 

F2  7 

-  *?F,  * 

(y+DbJDj  + 

2BoB;02  ♦  6B0 

FM  +  8B,F2S  +  4D,D12 

+  4BoB,D,2 

+  4D,D,n 
2  20 

+  (2y)bJD20  +  (2y)b'B6  + 

(^Kb, 

♦  (y+3)b’d 

1  9 

+  B  0 

2  13 

*  (2y)B„b; 

+  b(vi)b6d12 

+b1d,,  +  3 

(y-1)B>„ 

(B-117) 

+  (y-1)B 

,D,3  +  4<V- 

’  *B1D52  'Jh 

(y-1 )B  B  D 

VT  '  1  2  12 

F28 

=  8BiF,6  - 

«,8.D2  - 

- 

5B.D,,  +  (t-1)B,D„ 

(B-l 18) 

F29 

=  8B1F18  ' 

21^*216,8, 

Bs  -  4Bi01,  + 

(  y-1  )e,  b0  b, 

+  (y-UBjD. 

-  B  D 

3  5  Plul 

+  2  ^DSiB6 

Fso  =  (3y-1  )Bt  B2  83  +  4B2DJ6  +  (2y)B0B^  +  4B3D 


31 


Y+l 

2 


B1B3 


(B-119) 

(B-120) 

(B-121) 


32 


12B2Fi3  +  4B!Fi7  +  8Dn  +  6B22[32  +  ( 5Y" ^ ) Bj  B2  +  4(Y+1)B1B2D 


i  l 


+  W-DB2d14  +  (y-DBjD, 


25 


(B-122) 


146 


■  +  6B2F„  *  2B,F„  *  80„D12  *  D2Dm  *  !^»14 


+  2lY*l)BFDn  *  2(v+3)B1  B2Di3  t  2bJb2D2  *  (4>)B*B2 
*  (y-i)b2o13  *  (y-1)B,D24 


(B- 123) 


■  4B2F2  *  2(y+1)B1BiD1  *  4B0B2D2  *  12B0F,5  ♦  8B,F2]  *  80,0,, 

*  8Wn  +  4M,y  *  <2'KD,6  *  2(y*6Kb2  *  I^B2Dm 

+  2(y+1)BiB2Du  +  2(y-1)B2Di9  +  (y-DB0D)4  +  (y-DB^  (B-124) 


+  2(y-1)BiD34  +  4(y-1)B6Du  +  4(y-1)B0B1B2  +  2(Y-1  )B1  B2  Bfc 


=  2  (^l)6lBiB2  -  26j  D,  j 


(B-125) 


=  2(y+4)B 


(B-126) 


=  (y+1)b"B3  +  (y+1)B]B2B6  +  4Bj  Dj 6  +  4B2D19  +  (y)B^B2  +  (y+3)B0B'B2 

+  4B3D,2  +  4B.Du  +  2(y-1)B0D1,  *  l  (y-1)B0D14  *  2(y-1)B1D„ 

+  (y-1)BjDS4  +  2(y-1)B6D1  +  2(y-1)B0B1B6  +  (y-DBjB^  (B-127) 

+  (2y)B?D0  +  (y-DbX  +  4BqF2  +  2D;  ♦  4D0D2  +  4BoBlDl  +  2B>2 

1  r  5lr  3.1-  2  —  ,r  F3 1 

6  26  4  2  7  ^0  3  28^0  29  3  0  y  ~  2 

°  yo 

+  F32  [1  yo  (lnK)2  -  1  yolnyo  +  l  yo]  +  F33  [1  £  -  T6  \ 

+  F34  J  lYl*o  -  lyo]  +  F35[lny0  ’  ]]  +  \  F36^^ 


+  i  F37y, 


(B-128) 


_ _ 


147 
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1  „2, 


l  32  jui 


F 

+  1  d 

1  26 

8 

3y-l’ 

!  b2d 

6 

)  121 

1  n4  r 


+2Vd, 


(B-140) 
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CD3  =  CD30  +  CD31  +  CD32  +  CD3  3  +  CD34  +  CD35  ‘  (B‘323) 

Having  determined  the  (u  .Vj),  (u2,v2),  and  (u  ,v3)  velocity  components 
as  well  as  the  discharge  coefficient  constants  C  ,  CD2 ,  and  CD3 ,  the  important 
flowfield  variables  may  then  be  evaluated.  The  series  expansions  for  these 
quantities  are  given  in  Chapter  III  and  are  included  here  for  completeness, 

u(z,y)  =  l  +  u  =  l+u,e  +  u  c2  +  u,e3  +  ...  (B-324) 
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u(z,y)  =  1  +  U  =  1+U1C  +  U2£  +  Uj  £  + 


1  /2 

v(z,y)  =  v  =  £  [v!e  +  w2^2  +  Vj £3  +  ...  (B-325) 

M*(z,y)  =  (u2+  v2]l/2  =  1  +  Uj  c  +  u2e2  +  [u3+  v2  e5  +  ...(B-326) 
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While  the  preceding  list  of  constants  and  functions  looks  quite 
formidable,  it  is  to  be  noted  that  for  a  given  configuration  the  constants 
need  be  evaluated  just  once.  Only  z  and  the  functions  of  y  must  be 
determined  for  each  point  in  the  flowfield. 
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APPENDIX  C.  SOLUTION  SUMMARY  FOR  THE  CONVENTIONAL 
AX  I  SYMMETRIC  CONFIGURATION 


For  the  special  case  of  an  axisymmetric  nozzle  without  a  centerbody, 
the  parameters  which  are  to  be  specified  are:  y,  the  ratio  of  specific 
heats;  R  ,  the  wall  radius  of  curvature  non-dimensional ized  with  respect 

C 

to  the  throat  radius;  N,  the  number  of  terms  in  the  series  solution  to  be 
employed;  and  n  the  parameter  in  the  expansion  variable, 

c  =  [Rc+n]  1  »  (c-l) 

With  these  quantities  given,  the  perturbation  velocity  components  and  the 
discharge  coefficient  constants  for  the  first  three  solution  orders  may  be 
evaluated  from, 
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D3 

276,480 

The  quantities  u,  v,  M*,  6,  M,  p/pQ ,  and  CD  can  then  be  determined  from 
Eqs.  (III-73)-(III-78)  and  (III-80),  respectively. 
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The  resulting  solution  is  applicable  to  ax i symmetric  nozzles  with  either 
circular,  parabolic,  elliptic,  or  hyperbolic  arc  wall  contours.  It  should 
also  be  noted  that  the  expansion  parameter  employed  here  is  :  the  one 
obtained  by  simplifying  the  definition  of  e  used  in  the  general  solution, 

Eq.  (111-22),  for  the  case  of  a  straight  inner  wall.  The  latter  procedure 
results  in  the  definition  given  in  Eq.  (Ill -82 )  which,  because  of  the 
presence  of  the  coefficient  2,  is  somewhat  less  meaningful  for  the  present 
case  than  the  one  used  in  this  appendix,  Eq.  (C-l).  For  general  use,  and 
particularly  for  nozzles  with  a  small  wall  radius  of  curvature,  the  value 
n  =  1  is  recommended  for  the  solution  given  above. 


APPENDIX  D.  TABLES  OF  DATA,  FIGS  IV.11-IV.15 

COMPARISON  OF  FULL-  AND  HALF-SECTION  WALL  PRESSURE 
MEASUREMENTS  FOR  AX  I  SYMMETRIC  NOZZLE,  FIG.  IV. 11 


FULL-SECTION  MODEL 


Q  =  350.3  kPa  =  50.80  psia  Tq  =  288.9  K  =  520.0  R  v  =  1.4  (air) 
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HALF-SECTION  MODEL 

0  =  402.8  kPa  =  58.42  psia  TQ  =  288.9  K  =  520. 0°R  y  =  1.4  (air) 


h  „  =  80.0  k Pa  =  11.60  psia  ReOJ  =  3.09  x  10 
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FLOWFIELD  PRESSURE  TAP  MEASUREMENTS  FOR  CONVENTIONAL 
AXISYMMETRIC  NOZZLE,  FIG.  IV. 12 
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FLOWFIELD  PRESSURE  TAP  MEASUREMENTS  FOR  ANNULAR  NOZZLE 
WITH  CENTERBODY  CENTER  OF  CURVATURE  ALONG  Z  =  0.5  PLANE,  FIG. 


INTERPOLATED  CONSTANT  MACH  NUMBER  CONTOURS  FOR  ANNULAR  NOZZLE 
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